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Abstract 

We extend and develop our previous work on the evolution of a network of 
cosmic strings. The new treatment is based on an analysis of the probability 
distribution of the end-to-end distance, or extension, of a randomly chosen 
segment of left-moving string of given length. The description involves three 
distinct length scales: ^, related to the overall string density, ^, the persistence 
length along the string, and describing the small-scale structure, which is 
an important feature of the numerical simulations that have been done of 
this problem. An evolution equation is derived describing how the distribu- 
tion develops in time due to the combined effects of the universal expansion, 
of intercommuting and loop formation, and of gravitational radiation. With 
plausible assumptions about the unknown parameters in the model, we con- 
firm the conclusions of our previous study, that if gravitational radiation and 
small-scale structure effects are neglected, the two dominant length scales 
both scale in proportion to the horizon size. When the extra effects are in- 
cluded, we find that while ^ and f grow, C, initially does not. Eventually, 
however, it does appear to scale, at a much lower level, due to the effects of 
gravitational back-reaction. 
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I. INTRODUCTION 



Cosmic strings are topological defects that may be formed at a phase transition very 
early in the history of the Universe Because they are stable, they may survive to a 

much later epoch and thus provide one of the few direct links between the physics of the 
very early Universe and recent cosmology. In particular, they may play an important role in 
generating large-scale structure in the Universe . Observational tests of the idea include 
limits on the gravitational radiation emitted by collapsing loops and oscillating strings , 
gravitational lensing of characteristic form |llO|-|l^ and predicted anisotropy in the cosmic 
microwave background . 

All of these limits depend on our understanding of the process of evolution of a network 
of cosmic strings. One very important question is whether this evolutionary process leads 
to a "scaling" regime, in which the characteristic length scales describing the string network 
increase in proportion to the horizon distance P^^TB| . There have been numerical studies of 
this problem by at least three different groups [p!^^ -p!9|] . There is a wide measure of agreement 
between these groups. All find evidence for scaling of the large-scale structure. There is 
some disagreement over the predictions concerning the small-scale structure, but all groups 
agree that there is a substantial amount of structure on scales much less than the scale of 
the long-string network. It has certainly become apparent that a simple description in terms 
of a single scale is inadequate. 

In earlier work we sought to develop an analytic approach to this problem, to 

complement the numerical studies. In this work, we identified two distinct length scales, ^, 
related to the overall density of long strings, and ^, the distance over which the strings are 
correlated in direction. It became apparent, however, that this formalism is inadequate to 
represent the small-scale structure observed in the simulations. Accordingly, we have ex- 
tended our treatment by incorporating a third length scale, (, which describes the structure 
on the smallest scales. 

This treatment, which we present here, is based on a probability analysis of the string 
configuration, in particular the probability distribution of the end-to-end distance, or exten- 
sion, r, of a randomly chosen segment of string of length /; more precisely, of left-moving 
string. As before, we use null (characteristic) coordinates on the string world-sheet. In 
flat space, the left- and right-movers are completely decoupled; the universal expansion 
introduces a weak coupling between them. 

We should mention that an alternative approach to obtaining the evolution equations 
has been proposed by Embacher, in which a path integral formalism is used to obtain the 
probability distribution for the network of strings in flat space . This promising approach 
appears to produce similar results to our own in the circumstances where comparisons are 
possible. 

In Section II, we review the evolution equations developed earlier and introduce the 
fundamental probability distribution p[r(/)] on which our treatment is based. We aim to 
derive an evolution equation for this quantity, including terms representing the effects of 
stretching, gravitational radiation, intercommuting of long strings and formation of loops. 

For all except the very smallest scales, it is reasonable to assume that the probability 
distribution is Gaussian, characterized by the variance = K{1), say. Our evolution equa- 
tion for the probability p effectively reduces to an equation for this function K. The scales 
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associated with the large- and small-scale structure on the strings, ^ and (, are defined in 
terms of K, so their evolution equations follow from that for K. 

In Section III we analyze the Gaussian Ansatz in more detail and derive values of various 
expectation values that will be required in the subsequent analysis. In particular, we need 
a number of joint and conditional expectation values. 

The basic equations describing the effects of stretching, intercommuting and loop forma- 
tion are reviewed in Section IV. One of the very important things that emerges from this 
discussion is that these effects cannot properly be treated in isolation; they interact with 
each other in very complex ways. 

Another important lesson is the crucial significance of the correlations that develop be- 
tween left- and right-movers. These are discussed in Section V. In our previous work |]20|j2T|| 
we showed that the stretching process generates such a correlation, and derived a relation 
between the mean value a = — p q (where p and q are unit vectors along the left- and right- 
moving strings) and the length scale ^. More recently, however, we have realised that this is 
not the only significant process involved. In particular, loop formation also introduces cor- 
relations between the p and q vectors, by preferentially eliminating nearly matching pairs. 
The discussion of this process in Section V is based on an analysis of the angular probability 
distribution of these vectors. 

A similar process also introduces angular correlations on somewhat larger scales. This 
is important in determining the rate of loop formation. It is treated in Section VI. Section 
VII is devoted to the estimation of the parameters appearing in the stretching term. Then 
in Section VIII we deal with the effect of back-reaction from gravitational radiation, which 
operates on a very different scale from the other effects but is ultimately of great importance 
in the long-term evolution of the system. 

Finally, in Section IX all these terms are brought together to yield the overall evolution 
equations for the three length scales. The equations involve four dimensionless functions of 
the scale ratios as well as several dimensionless constants. Before analysing the solutions 
of the evolution equations, we study the behaviour of these functions in various regions of 
parameter space and try to estimate the constants. 

We then show that so long as the effects of gravitational radiation are negligible, the 
system behaves essentially as predicted by our previous study ||2l| . If we start with all three 
length scales approximately equal (and somewhat smaller than the horizon), ^ and later ^ 
will start to grow and, under reasonable assumptions, will evolve to a scaling regime in which 
they are approximately equal and proportional to t. The role of the important parameter q 
in our previous study is here played by a ratio between two of the dimensionless functions; 
it is, however, no longer a constant. 

There is another important parameter, k, also somewhat similar to q (or, more precisely, 
g — 1) but defined in terms of the small-scale structure. The value of k is crucial in deter- 
mining the behaviour of the third length scale (. If k were large enough, a complete scaling 
regime could be reached, with ^, ^ and ( all of comparable magnitude. This is not the type 
of behaviour seen in the simulations. A far more likely scenario is that k is less than the 
critical value. Then, while ^ and ^ oc t, the third length scale ( does not grow rapidly, and 
the ratio (/^ becomes very small. Eventually, when it is of order TGn, the gravitational 
radiation effect becomes significant and ( starts to grow. Here F is the numerical factor 
describing the efficiency of gravitational radiation, estimated to be of order 10^ 



3 



We then tackle the important question of whether a new scahng regime is reached, as 
various authors have aheady suggested |29-0l[|, due to the effects of gravitational radiation. 



It turns out that the answer to this question depends crucially on the value of another of the 
dimensionless constants we have introduced, (7, which relates to the effect of gravitational 
radiation on small-scale structure. Gravitational back-reaction tends to smooth out the 
small-scale kinks and thus to make ( grow; C defines the rate at which it does so. The 
essential condition for complete scaling is that C exceeds a critical value. 

We also consider the question of stability and show that if scaling is achieved it should 
be stable. 

The conclusions are summarized in Section X. We discuss in particular a number of 
remaining open questions and prospects for future study. 

II. EVOLUTION EQUATIONS 

We shall consider only the case of a flat Robert son- Walker universe, with metric 

ds^ = dt" - R\t) dx" = R^{r)[dT^ - dyi% (2.1) 

where R is the Robert son- Walker scale factor, r = x° is the "conformal" time, and x = {x^} 
are comoving spatial coordinates. 

A. The null coordinates 



For completeness, we recall here some of the basic formalism described in KG |]20|. It is 
convenient to use null (characteristic) coordinates u, v on the world sheet of the string. (In 
flat space, these are the coordinates u = t + a, v = t — a, where a is the length along the 
string.) We denote partial derivatives with respect to these coordinates by subscripts: 

dx dx . , 

Xu = ^, x, = — . (2.2) 

The null condition is 

xl = xl = 0. (2.3) 
In these coordinates the Nambu-Goto action is 

I = —jj, J dudv R^{t)xu-Xv, (2.4) 

where r = x^{u,v), and where the dot implies a scalar product in the Minkowski metric. 

In terms of the space-time coordinates (r, x) the equations of motion of the string may 
be written 

■^uv ~ ~^o{XuX^ + X^-X^), ^2 g'^ 

where 
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d.R IdR dR ^ 



where H is the Hubble parameter. 

It is convenient to define the unit vectors 



P = r^' q=r^' (2-7) 
which satisfy the equations of motion 

= -/ioX°(q- pp-q), , . 

q„ = -/lox;;(p-qq-p). 

Let us now consider a left-moving segment of string, i.e., a segment bounded by two 
values of the null coordinate m, say u\ and u^. 

The total physical extension, or end-to-end distance, of this segment at a given conformal 
time. To say, is 

RTdJ^^^ 



Ml 

R du 



du 



dv 



(2.9) 

In the second term, we can change variable to v and obtain 



rU2 rvi 

Ttot = -R / du:s.u-R dvyiy = l{ri- r^), (2.10) 

Jui Jv2 

say. (Recall that if Ui < U2 then vi > V2.) The two terms in ( p^.lUD are the left-moving and 
right-moving extensions, respectively. From now on, we consider only the left-moving term, 
and write 

r = ri. (2.11) 

The total physical length of the left-moving segment (measured along the string) may 
be defined to be 

l = 2R dux'l{u,v{u,To)). (2.12) 

J ui 

Note that because the time coordinates at the two ends have been chosen equal, we neces- 
sarily have, by analogy with ( p.lOj ), 

rU2 rv\ 

l = 2R duxl = 2R dvxl. (2.13) 

Jui J '02 

In this sense, therefore, the lengths of left- and right-moving string are exactly equal. 

It is useful to note that that the extension r can also be written in terms of the unit 
vector p, defined by (p.7|) , in the form 

'''dyp(y), (2.14) 

where y is a coordinate along the string defined by 

dy = 2Rxldu. (2.15) 
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B. Probability distribution of extension 



Consider a large comoving volume V, and let L be the total length of (left-moving) string 
within V. It is convenient to introduce the characteristic inter-string distance ^ defined by 

e = ^. (2.16) 

If we introduce a discretization scale 6, then L must be thought of as made up of small 
segments, each of length 6, with N = L/S. (There will be an equal length of right-moving 
string, but it is convenient to concentrate on one or the other.) 

Let us choose a particular length scale / and consider the probability distribution for the 
end-to-end distance (or extension) r: p[r{l)]d^r is the probability that a randomly chosen 
segment of length / will have extension r within the small volume d^r. Note that in contrast 
to KC |]2D| and CKA |^ we are here using real, rather than comoving, lengths, i.e., in terms 



of our previous notation, 

r = Ra, l = Rs. (2.17) 

The number of possible starting points within V is A^, so the expected number of segments 
with length between / and / + dl and extension r within d^r is 

jjp[r{l)]dh. (2.18) 

Because the distribution is highly non-random, there will be many segments with very 
similar values of r arising from overlapping selections, especially where at the end-points the 
orientations happen to be similar. 

Strictly speaking, some choices of starting point will not yield segments entirely within 
V; they will extend beyond the boundary, but will be matched by a similar number of 
segments entering V having originated outside it. 

Our object is, first, to derive an equation for the rate of change of p[r(/)], of the form 



dp 
di 





str. \ / er.rad. \ / l.s.i 




(2.19) 



and then to determine the nature of its solution. The various terms on the right repre- 
sent respectively the effects of stretching (due to the universal expansion), of gravitational 
radiation (back- reaction), of long-string intercommuting and of loop production. 

The separation between the intercommuting and loop-production terms is to some extent 
arbitrary: there is an upper limit to the size of a "large loop" . However, rather than imposing 
a sharp cutoff, we shall aim to count only those loops that survive the reconnection process. 
The reconnection probability provides a natural cutoff at a scale determined by the string 
density. There is no need to consider loops larger than this separately from the long-string 
network; such loops do not "know" that they are closed. They are rather likely to suffer 
reconnection. Their formation and reconnection may simply be regarded as instances of 
long-string intercommuting. 
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In addition to the equation (|2.19|) for p, we shall also need an equation for the rate of 
change of L: 



dL 
'dt 





gr.rad. 




(2.20) 



There is no term representing long-string intercommuting, which has no effect on L. 

The last term in ( |2.2(J| ), and the last two in ( [^.19D are each a combination of a negative 
term representing the effect of removal or destruction of segments and a positive term rep- 
resenting the corresponding creation of new segments. It is sometimes easier to consider the 
change in the expected number of segments, ( p.l8| ), due to one of these processes. It must 
be remembered that they affect L as well as p, but we note that from the rate of change of 
Lp we can easily find those of L and p separately, by using the normalization condition. For 
example. 




loops 



l^diLpY)^ 



loops 



and of course 




1 fd{Lp)\ 

^ V / loops 




(2.21) 



(2.22) 



loops 



The equations we obtain, not surprisingly, turn out to be very complicated. It is unlikely 
that an exact solution can be found. Our aim, however, is to find approximate solutions 
valid in special regions of interest, and in particular to establish whether a scaling solution 
exists. If it does, p should tend asymptotically to a limiting form 




t 



OO. 



(2.23) 



C. The Gaussian Ansatz 



For all except the smallest values of /, it is reasonable to assume that p is a Gaussian, 



p[r[ 



Here K{1) is the mean square extension, 

K{1) =^2 



3/2 



exp 



3 



d rr p[r{ 



(2.24) 



(2.25) 



From the equation ( p.l9|) for dp/dt we can derive a corresponding equation for dK/dt. 
We aim to show that its solution approaches the scaling form 



t — >• CX). 



(2.26) 
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For very large values of /, / S> t, we expect the string to behave like a Brownian random 
walk, so that K becomes a linear function of /, 

K{l,t) r^2^{t)l, />t, (2.27) 

with persistence length ^ (x t, but for smaller values of /, the variation of K with / will 
be more rapid, approaching ii' ~ as / — 0. (Of course, for such very small values, the 
Gaussian approximation breaks down.) 

In the next section, we examine the Gaussian Ansatz in more detail and derive some of 
its consequences. In particular, we evaluate various expectation values that will be needed 
in the subsequent analysis. Then in the following sections, we discuss the various terms in 
( pA9| ) and 



III. IMPLICATIONS OF THE GAUSSIAN ANSATZ 



We shall assume that except for very small values of / the probability distribution of 
the extension is a Gaussian, ( p.24| ). An important function, introduced in KG [^], is the 
correlation function 



/(!/)= p(0)-p(l/), 



(3.1) 



where y is the path-length variable introduced in (|2.15| ). This function is clearly related to 
the variance of r. In fact, we have 



K{l)^v^ = 2 dy{l-y)f{y), 
Jo 



from which it also follows that 



and 



K'{1)=2 fdyfiy), 
Jo 



K"{1) = 2/(/), 
where the primes denote derivatives with respect to /. 



(3.2) 



(3.3) 



(3.4) 



A. An illustrative model 

We shall not make any specific assumption about the form of as a function of I. 
However, it is useful to have a specific model in mind, to indicate the kind of behaviour we 
might expect. 

The simplest model, described in KG, is based on the hypothesis that there is a single 
scale in the problem: 

/(/) = e-™ (3.5) 
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for some constant D. (Note, however, that because of the change from comoving to real 
lengths (see Eq. ( p.l7|) ), D differs from the corresponding quantity in KC by a factor R.) 
Then 

i^(0 = -^(e-™-l+D/). (3.6) 
The correlation length ^ is defined by 

e= / dyf{y), (3.7) 
so that for large /, K is given by (|2.27|) . In the case of the single-scale model, 

e"=^- (3.8) 

This model does not describe very well the structure seen in the simulations. A better 
choice might be a model described by two scales: 

/(/) = (l-w;)e-^' + w;e-^', (3.9) 

where w is a constant in the range < w < 1 and we assume B. Here the small and 
large scales are 1/A and 1/-B, respectively. The expression for K is similar to (|3.6|) : 

m = ^^^(e-^' -l + Al) + |^(e-^' -1 + Bl). (3.10) 
In this case, the correlation length, defined by (|3.7|), is 

e = ^ + ---. (3.11) 

^ A B B ^ ^ 

Note that ^ is dominated by the large scale; the existence of the small-scale structure does 
affect it, via the constant w, but the size of the small scale 1/A is more or less immaterial. 

It is useful to define another characteristic length scale, (, related to the small-scale 
structure, by 



1 df{l) 



C dl 

For the single-scale model, C = but in the two-scale case 

1 



(3.12) 

1=0 



so 



^ {l-w)A + wB, (3.13) 



For the two-scale model, we may distinguish three distinct regions, in which the approx- 
imate forms of K (including first-order corrections) are 
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(^^^) R ^ '-^ - + 



ID K^wP-^ + ^-^- (3.15) 



Leading terms for very small and very large I may also be written in the forms 

(m) K^m---^. ^'-''^ 

We emphasize again, however, that this model is introduced for purely illustrative pur- 
poses; we make no specific assumptions at this stage about the form of K. 



B. Higher moments 

The requirement that the probability distribution is Gaussian means that all its moments 
are expressible in terms of the single function K. In particular, the variance of r^, the 
function 

K^2)il) =V¥ - (3.17) 

is given by 

K^,){l) = lK{ir. (3.18) 
The higher cumulants may be found from the cumulant generating function 

n=l ^• 

= -|ln(l - |2ir(/)). (3.19) 



C. Joint probabilities 

In evaluating the various terms in dp/dt, we shall encounter not only the probability 
distribution p[r(/)] but also various more complicated joint probabilities. We need therefore 
to extend our assumptions to cover these. 

Consider for example the configuration of two contiguous segments of lengths li and I2 
illustrated in Fig. 

We denote the joint probability of extensions ri and r2 within small volumes d^ri and 
(Pr2 by 

p[ri{h),r2{l2)]dhi£v2. (3.20) 

For a Brownian process, this probability factorizes into the product p[ri(/i)].p[r2(/2)], 
but this is possible only if K{1) is a linear function of /, since it implies that = + r"^. 
In the general case, where this is not true, TyT^ is non-zero. In fact. 
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K{h, h) = vT^= \[K{1) - K{h) - Kik)], (3.21) 

where of course 

l = h + l2. (3.22) 

It is interesting to examine the hmiting forms of the expression (|3.21|) in our two-scale 
model. The leading terms in the three regions are: 

it) /i,/2<i: K{h,h)^hl2] 

ill) i</i,/2<-I: K{h,l2) ^whh; (3.23) 
(m) /i,/2>:|: K{h,k)^^. 

The constant value as h and I2 approach infinity is noteworthy. 

It is useful to note that a similar formula to (|3.21|) holds for the mean value of the scalar 
product of overlapping segments; for example, in the configuration of Fig. |I|, 

rpr = |[ir(0 + K{h) - Kih)]. (3.24) 

This expression of course is non-zero even for linear K. 

In line with the Gaussian Ansatz for p[r(Z)], we shall assume that the joint probability 
( p.20| ) is also Gaussian, except of course for very small values of li or I2. Its form is then 
completely determined by the covariance matrix 



K = 



rrr2 

The distribution may be written 



where 



rrr2 

1-2 



K{h) K{h,l2) 



(det K)3/2 



exp[ 



^RK^R], 



R 



ri 
r2 



(3.25) 



(3.26) 



(3.27) 



The determinant in (|3.26|) may be written, using (|3.21|) , as a symmetric function of the three 
variables 



K = K{1), K,=K{h), K2 = K{l2), (3.28) 

namely 

det K = \[2KKi + 2KK2 + 2K1K2 - - Kj - K^]. (3.29) 

We note that the determinant is automatically positive provided that K increases as a 
function of / faster than linearly but less than quadratically, so that 

\2 



Ki + K2<K < {^Ki + ^K2Y. (3.30) 

This is assured if /(/) is a positive, monotonically decreasing function. 

It is useful to note that the joint probability ( 3.20 ) could equally well be written in terms 
of the variables ri and r, say, rather than ri and T21 provided of course that K were replaced 
by the appropriate covariance matrix. 



D. Conditional expectation values 



In later sections, we shall need various conditional expectation values, for example the 
expectation value of the extension ri of the li segment in Fig. ^ for a given value of the 
overall extension r of the composite segment of length /. Let us denote this conditional 
averaging, over the ensemble of segments with given values of r and /, by angle brackets. 

From the joint probability ( 3.26 ), it is straightforward to evaluate this conditional ex- 
pectation value, 

(n) = J d^nnpfnlr], p[r^\r] = (3.32) 

We find 



(ri) = (3.33) 
Kil) + Kih)-Kik) 

= — m) — ^ 

This is not at all surprising: a selected portion of a segment of extension r is obviously more 
likely to have its own extension ri in the same direction as r. 

For large values of / and h, for which K is linear, this yields, as one might expect 

(ri)f^yr (large /i,/). (3.35) 



From ( p.33| ) it follows of course that 



(r-n) = (3.36) 

Less obvious, but also useful, is the conditional expectation value of r^. A straightforward 
calculation yields 

(r;) = ^^--j^>% ^pp!r'. (3,37) 
From (|3.34| ) and (|3.37|) we obtain the conditional variance of ri: 

(r?) - ((n))^ = (3.38) 

It is remarkable that this conditional variance is in fact independent of the actual value of 
r. Note however that for any given value of r, it is less than the unconditional variance rf. 

A useful corollary of ( |3.34| ) yields the conditional expectation value of the unit vector p. 
We can write 

P(/i) = f^, (3.39) 

whence 
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v^{h) = \[K'{h) + K'{l-h)l (3.40) 

and 



(P('0)^^.^ "'"-'^f-'-' .. (3-41, 

For large values of / and /i, this reduces to 

(p(/i))^y (large /i,/), (3.42) 

again as one might expect. 

It is useful to note that these results, in particular ( |3.41| ), continue to hold for negative 
values of /i, provided that for negative /, K{1) is interpreted as meaning -ft'd/]). In other 
words, because the correlation of direction extends over a finite distance, a small segment 
close to but outside our chosen segment will still be correlated with it, though clearly less 
strongly than if it were inside. Of course, as li — oo, K'{li) —2^, while K'{l — li) 2^, 
so eventually (p(/i)) does approach zero. 



E. Triple joint probabilities 

We shall also need to consider more complex joint probabilities, such as the triple joint 
probability of the configuration shown in Fig. ^. 
This may be written as 



p[ro,ri,rj = ( — 



9/2 



(det K)3/2 



exp[-|R-K-iR] 



(3.43) 



where now 





"ro" 




R = 


ri 


, K = 




r 





ro-ri ro-r 



ro-ri 



ri-r 



ro-r ri-r 



(3.44) 



Integrating over all values of rg yields a joint probability distribution for ri and r identical 
in form to the previous one, save for the fact that the expression for rpr is different, namely 



rvr = \[K{lo + h) + K{1 - k) - K{lo) - K{1 - k - h)]. 



(3.45) 



With this change, the previous equations remain valid. 

Another useful result can be obtained from (|3.43 ). We can find an expression for the 
conditional average of a scalar product of unit vectors, (po-P2), where po = p(^o) and 
P2 = p(/o + h), by using analogues of ( p.39| ), namely 



Po 



P2 



dr2 



(3.46) 



with I2 = I — Iq — h- We find 
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(ro-rs) = ro-r2 + — pj^(r - H) 

= l[K{l) + K{1 -k- h) - K{1 - /o) - K{1 - h)] 

+ \m) + K{1,) - K{1 - l,)]m) + K{k) - K{1 - h)]^^. (3.47) 
Hence, differentiating with respect to both Iq and I2, we get 

+\[K'{lo) + K'{1 - mK'ih) + K\l - kf^^^. (3.48) 

Comparing with ( |3.41| ), we see that the conditional covariance function of po and p2 is again 
independent of r (but smaller than the unconditional value): 

(P0-P2)-(P0)-(P2) = |i^"(/l) 

- 4^[^'(^o) + K'{1 - loEK'ih) + K\l - h)]. (3.49) 
Note that as before these equations continue to hold even when either Iq or I2 is negative. 



F. Small-Z behaviour 

The Gaussian approximation obviously breaks down for very small values of I. To ob- 
tain information about the time evolution of the smallest-scale structures, we need to have 
approximate formulae that can be used in that region too. 

For very small /, the expectation value of takes the form given in (|3.16|(z)). 



^ - ^. (3.50) 

We shall also need to consider the probability distribution for a small segment of length 
li within a larger segment of length / and extension r. From (|3.41|) , it follows that the 



direction of ri is correlated with that of r. However, the length is not; the expression (p.37|) 



does not hold in the limit of small li. In fact, in that limit the length is essentially fixed; 
the probability distribution for ri is concentrated in a thin shell near |ri| = li. Therefore, 

{vD^vl^ll (3.51) 

Another important difference from the Gaussian case concerns the variance of r^, which 
is no longer given by ( p.l8|) . In fact, for small /, the leading term in K(^2) is clearly of order 
l^. To obtain a more specific result, consider for example a model in which small-angle kinks 
are randomly distributed on the string, with D kinks per unit length, and suppose the kink 
angles are distributed according to a (two-dimensional) distribution with (small) variance 
6'^. Then for values of / such that Dl <^ 1, we find 

K{1) =1^- \DI^{1 - ^ /2 _ ^niW. (3.52) 
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The characteristic scale ( of the small-scale structure, defined by ( p.l2| ), is 

2 

Wd 



C ^ (3.53) 



Similarly, we obtain 



K(^2){1) ^ ^Dl'>e\ 1^0. (3.54) 



30 

The ratio K(2)/K'^., which is | in the Gaussian case, becomes 

^(2)(0 ^ I 



K{lf 3002 ^ 

In particular, if the distribution of kink angles is Gaussian, then 

K{iy ^ 15 c' 



= -^Oas/^0. (3.55) 



(3.56) 



IV. DERIVATION OF BASIC RATE EQUATIONS 

The evolution of the system of strings is a complicated process. The mechanisms rep- 
resented by the various terms in ( p.l9| ) do not act independently. At least to a first ap- 
proximation, gravitational radiation is separable from the others, because the dominant 
scale involved is much smaller. We shall therefore postpone its discussion. As we shall see, 
however, the remaining three act in a complex synergy. 

We begin this section by deriving the basic equations for these processes, reviewing for 



completeness the discussion of KG and GKA pO| , ^ 



A. Rates of change of length and extension 



From the equations of motion ( |2.5| ), we can derive expressions for the expected rates of 
change of the length / and extension r of a chosen segment, expressing / and r as functions 
of / and r. Here, unlike KG, we use dots to denote derivatives with respect to real time; in 
particular, I stands for the average value, (dl/dt), of dl/dt over the ensemble of segments of 
given length and extension. 

Similarly, we can derive an expression for {L)stv., the contribution of stretching to the 
rate of increase in the overall length of string within a given large comoving volume. 

It is very important to note that there are consistency requirements. As before, let us 
denote by an overbar the average value of any function of r over the probability distribution 
p[r(/)] for given length I, e.g. 

1 = J dhip[r{l)]. (4.1) 

The lengths / of different segments will stretch by different amounts, depending on the value 
of r and other factors. However, the entire string can be chopped up conceptually into 
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segments of any prescribed length /, and its overall growth is obviously independent of /. 
Thus for every value of / we must have 



Deriving the expression for {dp/dt)stT. from / and r is actually quite subtle, because of 
the dependence of / on r; a sample of segments initially of equal lengths does not remain so. 
Consider a small time interval dt during which the expected changes of I and r are 

l^l' = 1 + idt, r ^ r = r + rdt. (4.3) 

The total length will also change of course, according to 

L^L' = L + {LUr.dt. (4.4) 

Now suppose that within V we select segments at random by choosing independently a 
random starting point and a random length /, chosen from a uniform distribution from up 
to a large upper limit, say L. The changes of length and extension of the chosen segments 
over a short time interval dt will vary randomly, with expectation values given by ([4.31). 



However, because of the consistency requirement ( [4. 2]) , the final distribution will still be 
uniformly distributed in /' (at least so long as / -C L, so that the upper cutoff of lengths is 
irrelevant). Hence we have the important equality 

p[r'{l'),t']^dh' = p[r{l),t]^dh. (4.5) 



Now 



and 



Hence it follows that 



dl' = ( 1 + dt^j dl (4.6) 



dh' = (l + dt^A dh. (4.7) 



This is consistent with the normalization condition because, from ([4.2|) , 

As before, let us consider a left-moving segment of string, defined by the inequalities 
Ml < u < U2- The length of our chosen segment at time t is 
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l{t) = 2R I duxyu,v{u,t)), 

ui 



while its extension is 



rU2 

r(t) = 2R / duy.u{u,v{u,t)) 

J Ml 

To find the rates of change, we use 



(4.10) 



(4.11) 



1 fdv' 

R [d^ 



dt 



Hence from the equations of motion ( pl^ 

dl R 



R [dv 



Rx?: 



(4.12) 



1 + 2 duxl 

J ni 



dt R 

R. 



R 



-1-2R rfMX^(l + p-q) 



Ml 



—2R I citix„p-q. 

•/ui 



(4.13) 



As before, we denote averaging over the ensemble of segments with given values of / and 
r by angle brackets. We then have 



fU2 

l = -2R du{x1p-q) = a{r,l)-l, 

Jui rC 



(4.14) 



say. Note that in order to satisfy the consistency requirement ( [4.2[ ), the average value of a 
over the r-distribution must be a constant a, independent of /. We must have 



(L)str. = 



(4.15) 



Thus a may be identified with the constant a of KG and CKA |^ . More generally, we 
shall find later that a has the form 



Similarly, we obtain 



a{r, I) = a + a{l) 



- K{1) 
K{1) ' 



(4.16) 



and hence 



ru2 /? 

r = -2i?/ dn(x°q)=/3(r,0-^r. 



(4.17) 



(4.18) 
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The parameters a and (5 are related but distinct. In particular, there is no special 
condition on the mean value of /?. It will emerge later that /? may be taken to be a function 
of / alone, independent of r. 

It is convenient to rewrite (|4.14| ) and ( [4.18| ) in terms of the path-length variable y intro- 
duced in ( p.l5|) . We then have 



al 



dy (p-q) 



(4.19) 



and 



/3r 



d?/(q). 



(4.20) 



We shall return to the computation of these averages in the next Section. 
Substituting / and r into ( [4.8|) , and using (|4.16| ), we thus find 



di, 



str. 



(4.21) 



where as before H is the Hubble parameter, H = R/R. 

Taking a moment of (|4.21| ), we find for the rate of change of K, 



dt 



-H 



d 



str. 



dl ( 



alK + al 



(2) 



K 



+ 2H(3K + HaK 



H 



2(3K - al 



dK d 



.K, 



dl 



dl K 



(2) 



(4.22) 



B. Intercommuting probability 



Next we review briefly the derivation of the intercommuting probability given in KC pO . 
This also provides an opportunity to refine the argument. 

Consider a large volume V, which, according to ( |2.16| ), contains a length L = Vj^ of 
long string. In other words, each volume ^ contains on average a length ^ of string. The 
model introduced in KC was to regard this string as formed of N independently moving 
straight segments, each of length where 

L V 

TV ^ - . -. (4.23) 

Choose any pair of these segments. For simplicity, assume that at the relevant time, 
the spatial coordinate a along the string is chosen to coincide with the variable y which 
measures length along the string. The two segments will intersect at some time during a 
short interval 5t if there is a solution to 

xoi + yix'i + t±i = xo2 + 2/2X2 + tX2, (4.24) 

with 
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0<yi<C, 0<y2<^, 0<t <St. (4.25) 

Equivalently, if the starting point of one of the segments, xqi, is fixed, intercommuting 
will occur if the other starting point, xo2, lies within a small volume 

6V = ^^5t|x; A X2-(ii - is)!- (4.26) 

(The additional factor of | appearing in KC, equation (4.26), was an error.) The probability 
of intercommuting between this pair of segments is 6V/V. To obtain the probability that a 
chosen segment undergoes intercommuting, we multiply by the number of segments, N, given 
by ( [4.23| ). Thus the probability that a string segment of length / will undergo intercommuting 
during a time interval dt is 

i dt / , »N 

X^, (4.27) 

where x is the average value of the scalar triple product \x.[ A X2-(xi — X2)|. 




The r.m.s. values of |x'| and |x| are y(1 + Ci)/2 and y(1 — «)/2, respectively. One way 
of estimating x is to use these as typical values. Then averaging over all angles, maintaining 
the orthogonality of x' and x, introduces a factor of 2/7r. Thus typically 

0.24. (4.28) 

Although X therefore has a weak a dependence, this is probably not sufficiently important 
to make it necessary to use an r-dependent value. 

Another way of estimating x is to rewrite the scalar triple product in terms of p and q 
vectors, as 

ilPi A P2-(qi - q2) + (pi - P2)-qi a q2|. (4.29) 

Averaging over all directions of the p and q vectors, assuming that they are independently 
and isotropically distributed, yields 

277- 

X = — ~ 0.18. (4.30) 



(It is easy to see why this result should be smaller than our first estimate (|4.28|). Using 



typical average values of the magnitudes of x' and x ignores the anticorrelation between 
them; including it would tend to reduce the estimate.) We can improve the estimate (|4.30|) 



by allowing for the correlation between the p and q vectors, corresponding to the non-zero 
value of a. This gives a corrected value 

X = - + ^ - 0.20. (4.31) 
^ 35 105 ^ ' 

This is probably the best estimate we can achieve. 

It might be argued that the model used here, assuming straight segments of length ^, 
is inaccurate. The length scale on which strings are roughly straight is not C, but ^. A 
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better approximation might be to assume that the string is composed of straight segments 
of length ^. However, this actually makes no difference to the final answer. The small 
volume 6V is then proportional to while N becomes N = V^/^^^, so the probability of 
intercommuting is still proportional to the length of the segment, in this case ^. 

The assumption that the segments are uncorrelated may seem rather dubious, particu- 
larly in the case where ^ is small compared to ^. Within each volume there would be 
a length C, of string, consisting not of a single straight section but of several segments of 
length ^. To treat these as moving independently is clearly to underestimate the number 
of intercommutings. However, the underestimate refers to intercommutings between nearby 
pairs, which will of course lead to loop formation and are therefore included elsewhere. 

One might also be concerned about the effects of small-scale structure on the strings. 
Segments that are roughly straight but kinky on a small scale would seem to offer a smaller 
cross-section. But in fact, this too is irrelevant to the probability of long-string intercom- 
muting, though very relevant to loop formation. Consider two kinky segments of string, 
each of length ^. For simplicity, suppose that one of them is instantaneously at rest, while 
the other is moving. During a short time interval 6t, the moving segment will trace out a 
thin ribbon of width |x|(5t. Intercommuting will occur if the other segment intersects the 
ribbon. So long as 6t is small compared to the small-scale structure, the probability of this 
is clearly proportional to the lengths of the segments and is in no way reduced by the fact 
that they are kinky. 

What is true of kinky strings is that where one intersection occurs, others are likely in 
the vicinity. This means that in the neighbourhood of an intercommuting loop formation is 
likely to occur. 



C. Effect of intercommuting 

We now turn to the effect of long-string intercommuting on the probability distribution 
p[r(0]. 

Clearly there is a negative term in {d{Lp) / dt)\ ^ ^ due to the destruction of segments by 
intercommuting, equal to 

- X^^Lp[v{l)]. (4.32) 

In this case, the corresponding positive term is rather more complicated. For each 
segment destroyed, a new segment of equal length / is created. However, its extension r is 
the sum of the extensions ri and V2 of the two, generally speaking uncorrelated, segments 
thus brought together. 

The number of segments created is exactly equal to the number destroyed. The proba- 
bility that one of these, chosen at random, is composed of lengths li and I — h, within an 
interval dli, is dli/l. Then the probability that the corresponding extensions are ri and r2, 
within intervals d^ri and d^r2 (see Fig. |^), is 

p[ri{h)]dhip[r2{l-h)]dh2. (4.33) 

To find the probability that the total extension is r, we have to set r2 = r — ri and integrate 
over ri. Thus the positive term in {d{Lp)/dt)i s,i. is 
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+ ^ fj^ j dh,p[v,{h)]p[{v-v,){l-h)]. (4.34) 

In contrast to the case of loop formation, here the assumption that the two probabihties are 
independent should be a good one, because the two segments involved generally come from 
regions that are far apart along the string. 
In this case, we evidently have 

(§),,-■ (-^' 

N / l.s.i. 

Hence, putting ( |4.32| ) and ( [4.34[ ) together, we obtain 
( dp\ X 

l.s.i 



/ « 1 3 



f^dh J dh,p[r^{h)]p[{r-r,){l-h)]~lp[r{l)]\. (4.36) 



It is again straightforward to find an expression for the rate of change of K. When we 
multiply by and integrate, in the first term it is best to go back to using ri and r2 as 
integration variables. Since there is assumed to be no correlation the mean value of ri-r2 
vanishes. Thus we find 



'dK\ X 

l.s.i. 



dt , . 

/ .S.I. ^ 



1^ dh[K{h) + K{1 - h)] - IK{1)\ . (4.37) 



Note that in general we expect 

K{1)> K{h) + K{l-h), (4.38) 

so the effect of long-string intercommuting is to reduce the value of K (unless i^' is a linear 
function of /, in which case it is unchanged). 

By symmetry, we can simplify the expression ( 4.37 ) slightly: 

fdK\ X 



dt , . e 

/ l.s.i ^ 



2 J^dhK{k) -IK{1)\ . (4.39) 



From this, we can find the rates of change of the various length scales. Intercommuting 
has no effect on L, or, therefore, on ^. 

The length scale ^ is defined by ( p.27| ). To find its rate of change we differentiate ( 4.39 ) 
and then allow / to approach infinity: 

2^ = „„, (4.40) 

dt l^oo dt ^ ' 



This yields 

Similarly, from ( |3.50|) we find 



(4.41) 
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(4.42) 



Of course, for this to be consistent, we also have to verify that the rates of change of K, K' 
and K" all vanish in the limit. This is easy to do. Thus we find 



Ci.s.i. ^ _xC 

c i 



(4.43) 



D. Effect of loop formation 

Next we review and revise the derivation of the probability of loop formation. This is 
perhaps the most problematic aspect of our analysis; we shall approach it in several stages. 

Consider a segment of string of length / and extension r. Wc want to evaluate the 
probability that this segment will form a loop during a short time interval. To be specific, 
let Q {r , I) dl dy dt be the probability that a loop will form in the time interval d,t, of length 
between / and / + dl, and with starting point within the small interval dy. This is the 
function we aim to estimate. 

From 0, we can compute rates of change due to loop formation of all the quantities 
we need. In particular, the probability that any particular point on the string will be 
incorporated into a loop within the time interval dt is clearly Xdt, where 

A= / dllA{l), (4.44) 
Jo 

and 

A{1) = J dhp[r{l)]Q{r,l). (4.45) 

The parameter A determines the overall rate of loss of length to loop formation. Consider 
a very long section of string, of total length L. Then the expected rate of change of L is 



— = -XL. (4.46) 




The dependence of © on r is unknown and may well be complicated. However, consistent 
with our Gaussian approximation, it seems reasonable to assume that, except for very small 
values of /, this is also Gaussian, parametrized by a variance function Q{1). Specifically, we 
assume that 

e(r,/) ^ A(/) (^:^Ii(^tMy'^'e-3rV2Q(0. (4.47) 

We now turn to the expression for the effect of loop formation on the evolution of ]9[r(/)]. 
Let us consider the expression for {d{Lp)/dt)ioops- It consists of a negative term representing 
the number of segments lost to loop formation and a corresponding positive one representing 
the number of new segments created by the process. 
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The chosen segment of length / and extension r will disappear if a loop is formed within 
it, overlapping either end, or enclosing it entirely. Let us ask for the probability that this 
happens within a short time interval dt due to the formation of a loop of length /i, within 
a range dli. 

We denote the starting point of the loop segment relative to that of our chosen segment 
of length / by uq (see Fig. ^). Allowing for all possible overlaps, the range of possible starting 
points is 

-li<y^<l. (4.48) 

[The total number of starting points, with discretization scale 5, is {l+li)/5] The probability 
that a segment of length li starting from one of these points has extension ri within a volume 
5^ is p[ri(/i)|r(/), ?/o]^^; where the symbol p[. . . | . . .] denotes a probability conditional on 
both the value of v{l) and on the position of the starting point Uq. Rather than using 
this conditional probability, it is more convenient to work with the corresponding joint 
probability, using the identity 

p[ri(/i),r(/)|i/o] =p[ri(/i)|r(/);i/o]p[r(0]. (4-49) 

However, the joint probability is still conditional on the position uq of the starting point. 

The probability that within a time interval dt a loop is formed of length between li 
and li + dli, starting between uq and ?/o + dy^ is, according to the discussion of Section V, 
9(ri, li) dli dyo dt. 

Thus we find for the negative term the expression 

-Lj^ dh j dhiQ{viM) j_^dyop[vi{h),r{l)\y^]. (4.50) 

In evaluating the corresponding positive term, we have to consider excision of a loop 
entirely within a segment of length I + Thus it is 

L dh dhiOinJi) dyop[riih),ir + ri)il + h)\yo]. (4.51) 
Jo J Jo 

Putting ([4.50|) and ( |4.51D together, we have 

(^^) ^^r^^^' I d'ri&{ri,h)yyyop[ri{h),{r + ri){l + h)\yo] 

- J^^^dyop[Ti{h),r{l)\yo]y (4.52) 
It is straightforward to perform the integration over r to find {dL/dt)\oops, because clearly 

Jdhp[ri{h),ril)\yo]=p[riih)], (4.53) 

independent of yQ. This of course reproduces ( [4.46| ) with A given by ( [4.44| ). 
Combining ( [4.52| ) and ( |4.46|) we obtain 

(^) " r^^^' Jdhie{ri,h)yyyop[ri{h),{r + ri){l + h)\yo] 

- J^^^dyop[ri{h),r{l)\yo] + hp[ri{hMm]. (4.54) 
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E. Probability of loop formation 



Now let us turn to the calculation of 9(r, I), or equivalently A(/) and Q{1). Our segment 
of length I will form a loop during a short time interval dt if the corresponding total ex- 
tension Ftot, given by (|2.10|) , vanishes at some instant during that interval. In other words, 
the corresponding left- and right-moving segments, each of length /, must have the same 
extension, r. 

To simplify the counting, let us imagine that space-time is partitioned into cells each of 
volume 5'^, and moreover that the string is partitioned into small segments of length 6. 

The probability that a left-moving segment of length / has an extension r within a 6^ 
volume labelled j is 

pj = p[v{l)]5\ (4.55) 

The probability that a loop is formed is essentially the probability that the corresponding 
right-moving segment also has the same extension. Assuming for the moment that the 
probabilities are uncorrelated (which as we shall see is by no means true), this probability is 
also pj. We have to multiply it by factors relating to dl, dt, etc. The number of length steps 
within a given range dl is dl/6. The number of time steps within dt is 2dt/6. (The factor of 
2 arises, as explained in KG pO| , because the segments are moving with the speed of light in 
opposite directions, so that each encounters a new segment after a time S/2.) The number 
of starting points on dy is dy/6. 

However, we should not merely multiply these factors. Particularly for small loops, the 
angles between the various vectors are small. 

A typical configuration is shown in Fig. |. The question is, by how much can we vary the 
length I, the starting point y and the time t without moving the difference in extensions of 
the two segments out of the 6^ volume? This is essentially the same calculation that we did 
earlier in estimating the parameter x that determines the rate of long-string intercommuting. 
The volume swept out by the total extension rtot when the parameters vary over ranges dl, 
dy and dt is 

|x'i Ax2-(xi -X2)|(i/c/?/cit. (4.56) 
Thus we should include in our expression for the loop-formation probability a factor 

A(r, /) = i(|pi A p2-(qi - qs) + (pi - P2)-qi A qsl). (4.57) 
The number of cells within the volume d^r is d^r/6^. Hence our required probability is 

p[r{l)]e{r,l)d'Tdldydt = p'^^j^^Air,l), (4.58) 

i.e., 

e(r,/) = 2p[r(/)]A(r,/). (4.59) 
Note that the factors of S cancel, as they must. 
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Inclusion of the factor A eliminates the problem of multiple counting of loops noted 



in KC [^. This was avoided in CKA |^ by imposing a small-scale cutoff and treating 
the contribution of small loops separately. However, that procedure introduced additional 
problems associated with the choice of the cutoff. It is better to treat loop production within 
a unified framework. 

So long as the loops are reasonably large, it is reasonable to assume that the unit vectors 
involved are independently randomly distributed on the unit sphere. Then we can replace 
A by its average, A = x ^ 0.2. 

For very small loops we expect A to be small. Consider for example the model in which 
the strings are composed of straight segments joining randomly distributed kinks. For very 
short loops one would expect to find only a single kink on each of the left- and right-moving 
segments forming the loop. However, it is easy to see that in this case the scalar triple 
product (|4.56| ) vanishes identically, because the three vectors are coplanar. This is because 
a triangular loop is necessarily planar. The point is that the loop-formation condition in 
this situation is satisfied only on a set of configurations of measure zero. 

This is of course an idealized model, but even in a more realistic model we should expect 
that the A factor would be very small for loops with only a single pair of kinks. 

It follows from this argument that A should vanish rapidly as / — 0. In the simple model, 
the leading contribution would come from loops with at least three kinks all together. If the 
kinks are randomly distributed, with a separation of order (, then the number of kinks on a 
segment of length / is Poisson distributed, with mean ^ so for three or more kinks we 
expect a factor of at least (//C)^- 

Also for small loops, the r-dependence is likely to be important. For values of |r| close 
to /, the segments forming the loop must be nearly straight with a high degree of correlation 
between the p vectors at the two ends, and also between the p and q vectors. In that case 
the scalar triple product will be small. On the other hand smaller values of r imply large- 
angle kinks and correspondingly less correlation between the unit vectors. So we expect A 
to decrease with increasing |r|. 

Except for very small loops, where A is in any case small, is generally much less than 
P. Thus it would seem reasonable to assume, consistent with our earlier assumptions, that 
it has the form of a Gaussian, 

A{T,l)^Aoe~^^"/^\ (4.60) 

where a is a constant of order unity. The overall factor Aq approaches x large / and 
vanishes at least like (l/C)^ as / — 0. The Gaussian factor yields a contribution to 1/Q{1) 
of magnitude 2a /3P. 

The formula ( [4.59D is still not accurate, for several reasons, but particularly because 



it neglects all correlations between the left- and right-moving strings. The rate of change 
of L due to loop formation is to be computed using (|4.46| ) and (|4.44| ). As it stands, this 



expression would give problems at both ends of the range of integration over /. At the upper 
end, we have to consider the probability of reconnection; we aim to include only loops that 
do not reconnect. At the lower end, the integral would diverge, because of the neglect of a 
very significant angular correlation effect. 

We deal first with reconnection. For a loop of size /, the probability of reconnection within 
a short time interval dt is x^dt/^"^. The probability that the loop will survive reconnection 
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to a much later time is therefore 




expi-x dt'—-^]. (4.61) 



If we assume that over the relevant period ^ at least approximately scales, i.e., ^(t) oc t, 
then the integral yields 

exp (-^) . im 

(It is not necessary here to allow for the variation of / with time, due to gravitational 
radiation, which occurs over a very much longer time scale.) To allow for the probability of 
reconnection, we should replace ( |4.59| ) by 

e(r, /) = 2p[r(/)]e-^'*/«'A(r, /). (4.63) 

The exponential provides an effective upper cutoff in ( |4.44| ) at a scale of order ^^/t. Above 
that scale, loops are rather likely to reconnect; below it, they mostly survive. 

It should be noted that the assumption of approximate scaling is a very weak one. Even 
if ^ does not exactly scale, the dominant contribution to the integral (|4.61|) will come from 



close to the lower limit so ([4.62| ) will change at most by a factor of order unity. 



Even with the modifications described, the expression (|4.63|) is still not entirely correct, 
because it neglects any possible correlation between the left- and right-moving segments. 
This turns out to be the most important effect of all. These correlations form the subject 
of the next Section. 

There is a potentially important correction yet to be included. Some of the long-string 
intercommuting events create large loops. (Small loops tend to be created by a rather 
different process, very dependent on the small-scale structure on the strings; there is less 
likelihood of confusion.) Unless these loop-forming crossings can somehow be excluded, 
there is a possibility of double counting. It is not immediately clear how serious the problem 
is likely to be, although one might guess that it affects positive and negative contributions 
in a similar way, so that the overall effect may not be very large. 



V. CORRELATIONS BETWEEN LEFT- AND RIGHT-MOVERS 

Taking account of the correlations between left- and right-moving segments, we should 
really write 

e(r, 1)p[t{1)] = 2p[r{iy, T{l)]e~^^'/^' A{r, I), (5.1) 

where p[r{l)]r'{l)]d^r d^r' is the joint probability that the corresponding left- and right- 
moving segments have extensions r and r' respectively. 

One approach might be to extend the Gaussian Ansatz, representing the joint probability 
distribution by a six- dimensional Gaussian with an appropriate covariance function r ■ r'. 
However, this is not a good representation of p[r(/); r'(/)]. The covariance is actually quite 
small (and negative), suggesting that the correlation effect yields merely a small reduction 
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in the loop-formation probability. But this is quite false: the joint probability distribution 
is in fact very sharply reduced in the forward (r = r') direction. 

We begin by considering the correlation of the individual p and q vectors. There are two 



quite separate processes that generate such correlations. In CKA we considered only 
the effect of stretching, but in fact loop production also plays an important role, and indeed 
intercommuting cannot be ignored. 

Consider a particular p segment, p(mo), and the approaching q segments, in particular 
a segment ci{vq) which encounters it at time to. (See Fig. ^) We wish to estimate how the 
angular probability distribution of q relative to the direction of p changes as the vectors 
approach one another. It is again convenient to use the path-length variable y defined in 
( p.l5| ). Let us define $(y,z), so that the probability that p(l/o)'q(z/i) = z, within the range 
dz, is 

Hy,z)^, (5.2) 

where y = yo — yi- With this normalization, the initial condition for $ at large yo — yi is 

$(00,2) = 1, (5.3) 

representing a completely random distribution when the segments are far apart. 
Let us now seek to write down an evolution equation for $, of the form 



(9$ 

dy 




(5.4) 



loops 



We are assuming here that the process of establishing an angular correlation as the vectors 
p and q approach takes a relatively short time, compared to the time scales for evolution, 
so that $ may be regarded as a function only of z and of the path-length y between them, 
not explicitly of the time. This is reasonable because the angular correlation sets in only 
when the vectors are already quite close. For similar reasons, we are justified in neglecting 
the effect of gravitational radiation which is small on the scales of interest here. 

We also introduce the corresponding angular probability distribution for the p vector 
that encounters q at yi. the probability that p(yo)'P(l/i) = z within the range dz is 

^ dz , , 

■^{y,z)—. (5.5) 



2 



Let us recall that 



/I 
^ — z^iy,z) = fiy), (5.6) 

where / is the function defined in ( p.l|) . (When it is necessary to distinguish, we denote the 
mean value of z with respect to the distribution $ by z^ and that with respect to \E' by %.) 
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A. The exponential Ansatz 



To reduce the problem to manageable proportions, we make a simplifying assumption 
concerning the angular distribution function ^, analogous to the Gaussian Ansatz, namely 
that (except when y is very small) it takes the form of an exponential: 

^(l/,^) = -^e^ (5.7) 

where /c is a function of y. The relation between k and / is 

^* = /(y) = coth/c-^. (5.8) 

(As we shall see the Ansatz breaks down for very small values oi y.) 

It will be useful to examine the limiting cases of large and small values of y. First, when 
y is large, f and consequently also k <^1. In that case we have 

k k^ 

/Ri--— (y large). (5.9) 

In most cases, therefore, a linear approximation k ^ 3f will be adequate. Neglecting k^, we 
may write 

^1 + kz !i:il + 3fz (?/ large). (5.10) 

In the opposite limit, where y is small (but not small enough to render the Ansatz 
invalid), we have 1 — / -C 1 and hence /c 3> 1. In this case, the leading approximation is 

k Y~~y ^ ~ small). (5-11) 

Here we may write 

^ ^ 2A;e-^(^-^) {y small). (5.12) 

The distribution becomes concentrated near z = 1 within a range of order 1/k. 

It is interesting to note that in the intermediate region, both approximations are in fact 
reasonably good. For example for / = |, the leading large-y and small-y approximations 
give respectively A; = | and k = 2; the correct answer is k = 1.8. 

When the exponential Ansatz is valid, all the moments of the distribution are of course 
determined by the expectation value. In particular, 

2 coth k 2 
= 1 ^ + 

+ — + — {y large). (5.13) 

Indeed, so long as A; <^ 1 it is a good approximation to set {z'^)q, ~ |. 

The exponential Ansatz can also be applied to the q distribution in the modified form 

*(^'^' = ssr6^""- 

(Since the directions of p and q vectors are an^z-correlated, with this definition b is positive.) 
However, as we shall see, the Ansatz ceases to be a good approximation for very small values 
of y. When the approximation is valid, the evolution equation for $ effectively reduces to 
an equation for b or equivalently z^. 
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B. Equation for $ 



Now let us consider the ensemble of approaching q vectors, in particular those that at a 
given time t fall within a small interval 6yi at yi. 

In the time interval dt, these vectors either move closer to yo or else are eliminated by 
being incorporated into a loop. The expected distance by which they move closer is 

dy = ~{2 + \y)dt, (5.15) 

with y = yo — yi- Here the 2 is the normal velocity of approach and the extra term Xydt 
represents the expected loss of length between yi and yo due to loop formation. Integrating 
this relation, we find 

2/ = ^ (e^(*-*) - l) , (5.16) 

where to is the time at which the segments coincide. Note that the distance y becomes 
exponentially large for time differences larger than 1/A. Loop formation effectively shuffles 
the string segments on this time scale, providing a very effective long-distance cut-off. 

During the time interval dt, the small interval 6yi on the string effectively becomes a 
little shorter, because some members of the ensemble are eliminated by incorporation into 
loops. In fact, d5yi/dt = —X6yi, so Syi is replaced by 6yi — Xdt 6yi. The angular distribution 
of the remaining vectors changes from $(?/, z) to $(?/ + dy, z), with dy given by ( |5.15| ). To 
find the difference between these, we need to know the angular distribution, X{y, z) say, of 
the vectors that have been eliminated. Here X(|/o — z)^ is the probability that a vector 
q which is incorporated into a loop within a short time interval has p(yo)'Cl(2/i) = z within 
the range dz. If we can estimate X, then we can obtain a differential equation for $. In fact 
(ignoring stretching and intercommuting contributions), 

5yi^y,z) = {6yi-Xdt6yi)<l>{y + dy,z) + Xdt6yiX{y,z), (5.17) 

or equivalently 

{2 + Xy)(^] =-A$ + AX. (5.18) 

V " / loops 

Note that this is consistent with maintenance of the normalization condition J = ^ 
provided that X is also normalized. 

The key observation is that since loops are formed by matching segments of left- and 
right-moving string, the angular distributions of excised p and q segments should be iden- 
tical. Clearly, the probability of obtaining a matching pair must depend on the probability 
distributions of both p and q. It seems reasonable to assume that X is proportional to the 
product \E'$. Normalization then requires 

X{y,z)=N{y)^{y,z)^y,z), (5.19) 

where 
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Thus we find 



/I //r 
^—m{y,z)^{y,z). (5.20) 

(2 + Ay) — =-A$ + AA^^$. (5.21) 

\ " / loops 

So far we have ignored both intercommuting and stretching, both of which will also have 
an effect on the angular distribution $. Consider first the effect of intercommuting. The 
probability that an intercommuting event occurs between yi and y^ within the time interval 
dt is xydt/^^- If it does occur the relevant members of the ensemble of q vectors are deleted 
and replaced by new vectors drawn from an essentially independent random distribution. In 
other words, $(?/, z) on the left-hand side of ( |5.17| ) is replaced by 

(l-^p?i)+^p?i. (5.22) 

Thus the effect of intercommuting is described by 

Finally, let us consider the effect of stretching. The equation of motion (|2.8|) for q may 
be written 

^ = -i/(pi-qq-pi). (5.24) 

where pi denotes the vector p(?/i) = p('Ui,fo) (see Fig. Strictly speaking, the Hubble 
parameter if is a variable quantity here, but, since the correlations extend over distances 
which are small compared to the horizon distance, it should be a good approximation to 
treat it constant. 

We now have to average this result over the angular distribution of p(?/i). In principle, 
this vector is correlated with both p(?/o) and q(?/i). However, any contribution proportional 
to q(?/i) cancels out in ( [5.241 ). Thus it is reasonable to consider only its correlation with 
p(?/o), and to replace p(?/i) by ^vi,p(i/o). 

There is an exactly similar expression for the rate of change of p(2/o)) which yields an 
identical contribution. Taking account of both we thus find 

d r n 

-2Hz^,— [{l-z^)^\. (5.25) 

str. \ / 

Bringing all three contributions together, we may write the equation for $ as 





9$ „ d 



xy . 



{2 + Xy)— = -2Hz^— (l-;z^)<l> - ^(1 - $) - A<1> + AiV^<l>. (5.26) 
oy oz J t, 

When the exponential Ansatz is valid, all we need is the equation for the rate of change of 
z^, namely 

(2 + Xy)^ = 2Hz4l - (F)^] + - Xz^ + Xzx, (5.27) 

where zx is the mean with respect to the distribution X = A^\E'<1>. 
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C. Equation for ^ 



Before trying to solve the equation for $, some parenthetical remarks about the possi- 
bility of deriving a similar equation for \[' may be in order. 

The effects of stretching and intercommuting on \E' are very similar, and we have assumed 
that the distribution X of excised segments is the same for both. So we should be able 
to write down a very similar equation for the evolution of \E'. However, there is a very 
important difference, concerning the rate at which one p segment approaches another. The 
most obvious difference is that the distance y = yo — yi between the two p vectors decreases 
only because of loop formation. In other words the term 2 in the factor on the left of ( ^.18[ ) 
is absent in the corresponding equation for ^ . 

But there is a further point: it is no longer reasonable to assume that dy/dt is independent 
of z. In fact, if it were, the string would never develop the long-range directional correlation 
that it does. The formation of loops (except for the very smallest) depends strongly on the 
large-scale configuration of the strings. If the left- and right-moving sections are relatively 
straight, the large- loop formation probability is low, because the values of r for given / are 
large. On the other hand, if the strings are curled up tightly, typical values of r are small 
and large-loop formation becomes highly probable. In the former case, positive values of 
p(l/o)'P(l/i) are clearly favoured. Conversely, if p(yo)'P(2/i) is positive, the expected value of 
—dy/dt will be smaller than if it is negative. Roughly speaking, we may expect 

^^^-[\y-{z-z^)v{y)l (5.28) 

where v{y) is a function that could in principle be determined from the later discussion 
of the detailed effects of loop formation. There is no doubt a similar effect even for the 
cross-correlation between p and q, but in that case it seems likely to be negligibly small. 
It follows that for \E', in place of ( p.26| ), we would have 

— [[\y -{z- -zs,)v{y)\^) = -'^Hz^— [(1 - z^)^] - ||(1 - ^) + AiV$^. (5.29) 



Similarly, the analogue of ( |5.27|) is 



— (Ay^^ - [{z^)^ - ziHy)) = 2Hz4l - (z^)^] + ^z^ + A%. (5.30) 



D. Solution of equation for $ 

Having set up the equation (|5.26| ) for $ we now set about solving it. For the moment 
at least, we shall treat as given, via the exponential Ansatz, in terms of the correlation 
function = f{y). 

Consider first the region of large y where the linear approximation ( |5.10 ) for \E' and $ 
should be valid. In this case, we need consider only the evolution equation for So long 
as Zip remains small, we can also use the approximation (-2^)$ ~ ^- Note also that under 
these conditions, the distribution X is also linear, with 
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zx = Zil, + z<s>. (5.31) 

Thus we find 

(2 + \y)^ = + (a + ^ff) -z^. (5.32) 

It is straightforward to integrate this equation using an integrating factor, to obtain 



where 



n = ^. (5.34) 



As we noted earher, the exponential Ansatz breaks down near y = 0, but if we ignore 
that for the moment, we can estimate the value of a = —z<s>{0), as 



a ~ 



fdy[^ + Yl ^'""'^'^'^M- (5.35) 



For large values of y, z^, is expected to fall off exponentially. The model described in 
Section 2 suggests that z^^ ^ we'^^, with B = w/^. In that case, z$ is expressible in terms 
of the incomplete gamma function: 

z^(y) = - '^^^^^^^ x-"e"-^^r(n,x) (5.36) 



x=[^ + ^){2 + Xy). (5.37) 



where 



A 2 

Expanding for large x, the leading term is 

(X + ^H)we-^y 

As a consistency check, we may substitute this solution into the equation ( |5.3CI| ) for 
the laxge-y \& distribution and verify that it can be satisfied with a reasonable form of the 
unknown velocity-distortion function v{y). In fact, we find for the leading approximation 

v{y) = ^[w\ + ^^ye~''y, (5.39) 

which seems entirely reasonable. 

For y ^ 1/X, —zq> is small compared to z^,. But as y falls it grows rapidly. For moderate 
values of y the two are of comparable magnitude, assuming that 1/ A is of the same order as 
i and i. 
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We know from the simulations that even at y = 0, does not become large P,p!9|: 

^,^(0) = -a^ -0.14 (5.40) 

(in the radiation-dominated era). Hence it is reasonable to assume for all values of y that 
z$ <C 1. This does not necessarily mean that the exponential Ansatz is valid (a point we 
shall return to shortly), but so long as it is we can still use ( [5.27| ), with (z^)$ ~ |. However, 
as ?/ — > 0, we can no longer use ( 5.31 ); instead, we have both ^ 1 and — ^ 1- The net 
effect is that in the small-y region, the exponent n — 1 in ([5. 33]) becomes n. But since Xy 



and Xy' <^ 2 in that region, the effect is minimal. The expressions obtained above should 
still give a good approximation to 



E. Behaviour near y = 

To get at least a rough estimate of the value of z$(0) = —a, let us first assume that 
^ 1 and ^ 1. Then in ( ^.36| ) both x and n are small. The leading term in the 
expansion of T{n, x) for small x and n yields 



q; ~ w 



l + i|)ln(Ae) (Ae»l). (5.41) 



This is, however, much too large, clearly inconsistent with our assumption that % always 
remains small. 

It is perhaps more plausible to assume that A^ <^ 1 and A^ <^ 1. Then both x and n are 
large, and roughly equal, and the asymptotic form of the incomplete gamma function gives 

(A^ < 1). (5.42) 

Although z^ remains small, the exponential Ansatz is not in fact a good approximation 
near y = 0, because the distribution \l/ becomes so sharply peaked near z = 1. The distribu- 
tion functions in the last two terms of ( p.26| ) are of course both normalized, but the second 




one is negligibly small over most of the angular range, becoming very large near z = 1. 
Thus, while a linear approximation to $ remains good for most values of z, near z = 1 it 
becomes very poor. As y — 0, $ acquires a deep hole in the forward direction. For this 
reason, our estimate of ^$(0) requires some correction. 

To estimate the likely size of the effect, let us assume that the exponential Ansatz is at 
least qualitatively reasonable down to values of y of order (. So when we come to consider 
the equation ( |5.26 ) in the region of small y we can use the exponential (or indeed linear) 



form as our initial condition at ?/ ~ 

When y ( we can use the approximation ( |5.12|) for \1', with zq, ^ 1 — {y/(). Then our 
equation (|5.26| ) for $ becomes 



(2 + Xy)^ = -2H (l-A-^ [(1 - ^2)$] _ ||(i _ $) _ A$ + XN'^e'^^-'^'^^y^. (5.43) 

It is convenient to change from y to k ^ (/y a.s the independent variable. The equation 
involves three small parameters, X(, H( and xC^/^^- Since the last of these is very small 
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indeed, the intercommuting term will give a very small contribution in this region. Moreover, 
we expect HC, <C AC -C 1, so the contribution of the stretching term is also likely to be small. 
Neglecting these terms, and also neglecting \y in comparison with 2, we find 

= AC(1 - 2Arfce-'=(^-^))$. (5.44) 

OK 

We now seek to solve this equation, starting with some initial value feg of k, at which we 
assume that $ has the exponential form (|5.14|) , with a value ho of h. Since XC, is small, we 
may adopt a perturbation approach, taking as our zero-order solution, ^^{k.z) = ^{kQ,z), 
independent of k. Correspondingly, the zero-order value of N is given by 



7' = t ^'^ke-'^'-Mh. z) ^ $(fco, 1) = (5.45) 
J-i 2 smhoo 

Now we substitute this into the right hand side of (|5.44 ) and integrate, obtaining in first 
order 



*,.,.)^*(.„,.)e.pg-|-ACiV„£fe-..-.) 



= *(*:„, j)exp ('^ - ^ + \CN„[Ei(-h,{l - z)) - Ei(-<:(1 - z))]\ , (5.46) 

where Ei is the exponential integral. 

It is interesting to examine the special case z = 1. At that point, we have 

m 1) ^ N,' exp f ^ - ^) • (5-47) 



Note that $(/c, 1) approaches zero as /c — * oo (or y —>■ 0), but very slowly. 

The most interesting limit of course is A; ^ oo. Strictly speaking, in this limit the first- 
order approximation in A^ breaks down, but in fact it does so only at such large values of k 
that d^/dk is already negligibly small. There seems no need to go beyond first order. For 
2; 7^ 1, we have in that limit 

^k = 00, ^) = $(fco, z) exp + ACiVoEi(-A;o(l - z))j . (5.48) 

The behaviour of this function near 2; = 1 is given by 

^{k = oo,z)^ Nq^c-^"^^-'^ exp 1^^^ [ko{l - z)]^'^^°, (1 - ;z < 1). (5.49) 

This clearly exhibits the expected sharp dip near z = 1: the last factor vanishes at that 
point, but is close to unity over most of its range. 

We are particularly interested in the average value ^$(y = 0) = —a, obtained by mul- 
tiplying ( |5.48D by z and integrating. Since $(?/ = 0,2;) is a product of two factors, each of 
which separately gives a small value of z, the two effects are approximately additive, and 
we may write 
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a 



aiin + ttnl, 



(5.50) 



where the subscripts stand for 'hnear' and 'non-hnear'. Here the hnear contribution is 
the value given by ( |5.41| ), ( p. 421 ), or something in between, while the non-linear term cini 
comes from the second factor in ( |5.48|) . To first order in AC, it is 



2 ■ 



(5.51) 



VI. ANGULAR CORRELATION AND LOOP FORMATION 

We are now in a position to return to the calculation of the probability of loop production, 
Q{r,l). The provisional formula (|4.63|) is wrong because it assumes that the probability 
distributions of the left- and right-moving segments are independent. It should be replaced 
by (13). 

To complete our programme of expressing all the terms in the equation for dp[r{l)]/dt 
in terms of p itself, we have to reexpress the joint probability here in terms of individual 
probabilities. 



A. Small loops 

Let us first consider the case of very small loops. We shall deal separately with larger 
loops in the next subsection. 

In the case of small loops, the essential effect is due to the angular correlation between p 
and q vectors. To be completely correct, we should consider the joint probability distribution 
of all the p and q vectors forming this section of string. This would obviously be a very 
complicated object; we are not in a position to deal with it. However, in the case of small 
loops, the internal correlation between p vectors at different points is very strong (as is 
that between q vectors), so it seems reasonable, in order to represent the effect of the p-q 
correlation, to choose a single representative vector from each class. We choose the pair for 
which the effect is strongest, namely the vectors at the mid-points of the segments, namely 
p(//2) and q(//2). (We have considered as an alternative averaging over the chosen position; 
this makes little difference.) 

Thus we take 

p[r(0;r(/)] ^ |^^<|.(0,p-q)p[r(0|p(//2)]p[r(0|q(V2)], (6.1) 

where p[r(/)|p(//2)] is the probability distribution of extension conditional on the direction 
of the vector p at the mid-point, d^p denotes an integration over the unit sphere and of 
course $ is the angular distribution function for z = p-q evaluated at the point where the 
two vectors meet, namely y = 0. It is reasonable to assume that, apart from their mutual 
correlation, the p and q vectors are uniformly distributed on the sphere. 

The Gaussian approximation should be valid, both for the conditional probabilities and 
for the loop-production function 0(r,/), since there is an effective cutoff in the region of 
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very small loops due to the behaviour of the angular distribution function and the volume 
factor. Then p[r(/)|p(//2)] is the Gaussian distribution with appropriate values of r and r^. 
By ( CT) , 

r = ^p, (6.2) 

where 

S = rp = K'{l/2). (6.3) 

The value of is unaffected by p: = K{1). This means of course that the variance of r 
is reduced: 

K{1) -f^ = K - S^. (6.4) 

It is worth remarking that for moderately small values of / the Gaussian Ansatz is 
a much better approximation for the conditional probability p[r(/)|p(//2)] than it is for 
the unconditional probability j9[r(/)]. For small /, the probability distribution is of course 
concentrated near the sphere = P and is nothing like a Gaussian centred at r = 0. 
However, the Gaussian approximation to p[r(Z)|p(//2)] is centred near that sphere and has 
much smaller variance, since for small / 

S^l-^, (small/). (6.5) 

We can now compute G, or the variance function Q in the Gaussian approximation. 



Substituting (pTTp into (|5.1|) , we obtain 

3/2 



e(r, i)pm ^ Ail) (^g§|y|f^) e--^[^/-«-V2Qa)] 

^2e->^'*/«'A(r,/) | ^ ^ $(0, p-q) p[r(0 |p(//2)] p[r(/)|q(//2)]. 

(6.6) 

It is easy to carry out the integrations over p and q, leaving only a single integration which 
may be written 

e(r,/)p[r(/)] ^2e-^'*/«'A(r,/) (^_^^\-Hr^+s-)/K 

^. ; J..^fn 0..2 i\„.„t, / 



j dM$(0,2M2- l)sinh (^-^uj . (6.7) 



3S'r Jo 

Without having a specific form for $, it is not possible to proceed further. However, we can 
get a good idea of the likely effect by assuming that $ has a sharp step-function cutoff, i.e., 
$ = for z > zq = 1 — 2ani and $ = constant for z < zq. We then find 

e(r,/)Mr(/)] ^ 2e->^'*/^^A(r,0 (^)' e'^^^^^^^^^^^^, (6.8) 



36 



with 



3Sr 

k ■ 



(6.9) 



Using ( [4. 601) , we can now perform the integration over r in ([4.45|) to find A. For small /, 
S'^ / K 6(/l, so there is a strong exponential cutoff for / < (. Thus we can legitimately set 
A(r, /) ^ Aq. Then we obtain 



A 



Ar 



47r(l - ani) \ttKJ S 



1/2 1 



-3S'^/K 



(6.10) 



B. The intermediate-scale region 

The formula ( |6.1|) incorporates our model of the p-q correlation which, as we have seen, 
provides a very effective cutoff for values of I of order ( or less. However, it would still 
predict an impossibly large value of A. If we substitute from (|6.10| ) into ( [4.44| ), this arises 
from the contribution of the intermediate region, where 

c«/«e. (6.11) 

In this region, K ~ w{l — w)l'^, while S'^/K ^ w/{l — w) is of order unity, so the integrand 
in (|4.44|) behaves like Thus we find A ~ 1/C, up to a constant of order one. This 



would mean an extremely rapid decrease of L, on a time-scale of order which is clearly 
not consistent with the results of the simulations. 

The explanation for this discrepancy again lies in the angular correlation effect, but of 
longer segments of string, not merely individual p and q vectors. 

Think of a section of string of length / and extension r, and suppose that C ^ ^ ^ ^- In 
other words, the section contains many kinks, but viewed on a large scale it is likely to be 
fairly straight, i.e., |r| is a sizable fraction of /. Now consider the collection of p vectors on 
this section. Their ensemble average is of course (p) = r/l. So the distribution of p vectors 
will be strongly skewed, concentrated into a cone around the direction of r. 

The q vectors on the corresponding right-moving section will also be concentrated in a 
cone, around the direction of r^. In most cases, the two cones will not overlap much, and 
rather few loops will form. On the other hand, where the cones do overlap, many loops will 
form. Such a section will disappear rapidly in a burst of loop formation. After a short time, 
regions where the distributions overlap will be rare. 

Of course, regions of overlap are continually being reformed, as sections of string meet 
new partners. In particular, the intercommuting process brings together segments of string 
that had been far apart and are therefore more or less uncorrelated. In some fraction of 
cases, depending on the size of the cones, these segments will have overlapping distributions, 
so a fresh burst of loop formation will be triggered. This burst phenomenon has indeed been 
observed in the simulations. 

It is important to note that the angular concentration effect will be even more pronounced 
for somewhat larger loops. One may think of the string section as composed of short straight 
segments of length ~ C, with randomly varying orientation within the cone around r. Thus 
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their transverse extensions are essentially a two-dimensional random walk. If we select a 
length nC, its overall extension will be close to n^r/l, with a transverse spread proportional 
to \/n(. Hence the angular distribution of such sections will be concentrated in a cone whose 
angle is reduced by a factor 1 / y/n compared to the cone of p vectors. 

The question we want to answer is: how does the loop formation probability depend on 
the three scale lengths C;^?^? 

In approaching this question, let us begin by considering an initial state in which all 
three scales are comparable in magnitude, as would be expected shortly after the string- 
forming phase transition: C ~ ^ ~ ^ < ^- this case, the upper cutoff of the loop formation 
integral, at C.'^/xt^ may well be smaller than the lower cutoff, (. What this means is simply 
that almost all loops formed reconnect to the network. There is then no reason why either 
( OT should grow rapidly. If we consider only the stretching terms, we expect (in the 
radiation era) 

i = ^/f=^. (6.12) 

Here a may be dependent on the ratios of length scales, but its value is presumably small, 
say around 0.2. This implies that C, oc t^'"^. Thus we reach a regime where ( ^ ^ ^ < t. 
The upper cutoff does grow, if only rather slowly: C.'^/x't oc t^ "^. 

The next stage is easy to describe, at least qualitatively. Large numbers of small loops 
start to form, without reconnection, increasing the rate of growth of ^. But also ^ starts 
to grow. This is because of the selective nature of loop formation. Some sections of string 
will be relatively straight. On those, the distributions of p and q vectors will be confined to 
cones. In a few such cases, the cones may overlap, leading to rapid disappearance; in most, 
they will not. By contrast, on the more curled sections of string, p or q vectors will readily 
find partners. Those sections will disappear too, though not perhaps quite as rapidly as 
where there are overlapping cones. The net result will be to eliminate selectively the more 
curled sections of string, leading to growth of ^. 

One consequence of this that will be important later is that the function 0(r,/) rep- 
resenting the probability of loop formation will be more concentrated towards r = than 
would otherwise have been expected. In other words, the variance function Q defined in 
( [4.47| ) will be smaller than the original forms ( [4.59|) and (|4.6CI| ) would suggest. 



Although ^ starts to grow, there is still no reason for ( to do so. Every intercommuting 
event and every loop formation introduces new kinks, keeping down the average interkink 
distance. (As we shall see later, there is good reason to believe that ( does eventually start 
to grow, but only when gravitational radiation becomes important.) 

Thus, within a few expansion times, we may expect to reach a regime where C ^ ^ ^ < 
t. How, in such a regime, does the loop formation rate depend on the various length scales? 
In other words, what correction factor do we need to apply to the rate (|5.1| ) calculated 
neglecting the angular correlation effects? 

Consider a particular loop size / within the range ( '€i I ^"^/xt- We also assume 
that / -C ^. Typically then a segment of length / will be part of a longer section that is 
roughly straight on a scale of order ^, with many kinks separated by distances of order (. 
The segments of length / will be concentrated around this overall direction. Unless this is 
a region where there is a burst of loop formation, the corresponding right-moving segments 
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will be concentrated in a similar, but essentially non-overlapping, cone — simply because 
regions where cones overlap quickly disappear. 

In these circumstances, loop formation will occur only when segments meet new partners 
with a different overall direction. This happens for two reasons — intercommuting and the 
steady progression of left- and right-movers. 

Consider first intercommuting. The probability that a segment of length I experiences 
intercommuting within a time interval 6t is x^^'^/^'^i o^') to put it another way, the average 
time between such events is ^^/x^- 

The steady progression will bring a given segment alongside ones with a different overall 
orientation when the left- and right-movers have moved relatively by a distance of order ^, 
i.e., after a time ^/2. In the regime we are presently considering, this is a much shorter 
time, so the steady progression is the important effect. (The intercommuting time would be 
shorter only for segments well above the upper cutoff length.) It is therefore reasonable to 
ignore the intercommuting effect. 

It is important to realise that this does not mean intercommuting plays no role; quite the 
contrary. What we are saying is that the direct effect of intercommuting within the chosen 
segment is unimportant. The indirect effect of nearby intercommuting is one of the things 
that keeps ^ down, and is in fact clearly crucial. 

If the orientations of the segments of length I were random, a given segment would meet 
another of essentially different orientation every time it moved on by a distance /, i.e., after 
a time 1/2. As it is, because of the long-range correlation, it will meet a segment with 
an essentially different orientation only after a time ^/2. On this ground, the rate of loop 
formation should be suppressed by a factor of approximately l/^ (for I < C,', there is no 
suppression for / > ^). 

This is not the end of the story, however. We have to consider what happens when a new 
segment is encountered. The new right-moving section of string may have any orientation 
relative to the left-moving section. We may assume that the relative orientation is random. 
Thus if a is the angle between the two, then cos a should be uniformly distributed between 
— 1 and 1. 

As we saw, the left-moving segments of length / lie within a cone whose semi-vertical 
angle (3 is proportional to 1/ ^/n, where n is, roughly speaking, the number of sizable kinks 
on the segment, namely n ^ The solid angle within the cone is of order vr/?^ ~ ttC/^- 

Now clearly, if a is significantly larger than /3, then our chosen segment is very unlikely 
to meet a matching partner. On the other hand, in the relatively rare cases in which a < /3, 
loop formation is extremely probable. In other words, we may expect another suppression 
factor, roughly equal to the solid angle of the cone divided by Air, namely C/4/. 

Putting these two suppression factors together, we obtain an additional factor in A(Z) 
of about We now insert the expression for A(/), given by ( |6.10| ) with this factor, into 

the formula ( ^.44| ) for the overall loop-formation rate A. Here we have to integrate from the 



lower limit ( to the upper limit ^^/(x^)- Since A(Z) oc l/l^, the integral behaves like / dl/P, 
and it is the lower limit that dominates — as must be true if most loops produced are indeed 
small. Hence the integral is of order Taking account of the extra suppression factor, 
we see from ( |6.1CI| ) and ( [4.44| ) that 



A = |, (6.13) 
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where E is a dimensionless parameter rather less than unity. 
Note that the index n defined in ( 5.34|) is 



So when ^ <^ ^, n may be small, but for ^ ~ ^, it will be of order unity or even somewhat 
larger. 



C. Rates of change of ^ and ^ 



We now take up the calculation of the effect of loop formation on the rates of change of 
the length scales in the problem. 

First, we consider ^. It is of course defined in terms of L by ( p.l6| ), from which it follows 
that 



3^ 



L 

2L' 



(6.15) 



Thus the rate of change of ^ due to loop formation is 



Cloops 



A 
2 



E 
2?' 



(6.16) 



where the dimensionless function E is given by 

Ae = n dllAil). 
Jo 



(6.17) 



For the other two scales, we have to examine the rate of change of the variance function 
K. The Gaussian approximation should be valid, since there is an effective cutoff in the 
region of very small loops due to the behaviour of the angular distribution function and the 
volume factor. 

It is then straightforward to derive an expression for the rate of change of K by performing 
the integrations over ri and r in ( [4.54|) explicitly. 

In the Gaussian approximation, the joint probabilities are given by expressions of the 
form ( |3.26|) with K replaced by the covariance matrix ( |3.31|) and FrT given by ( p.45|) . 

For the first term of ( [4.54| ), the appropriate covariance matrix is 



K 



(6.18) 



where 



(6.19) 
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while for the second it is 



K 



K 



(6.20) 



where 



= \[K{h + yo) + K{1 - 2/0) - K{yo) - K{1 -1,- y^)], 
K = K{1). 



(6.21) 



The exponent in the first term of the integrand in ( [4.54| ) contains the inverse of (|6.18|) 
plus a single-entry contribution from the extra factor G(ri,/i). In other words, the inverse 
of K is replaced by the inverse of a new matrix L, given by 



-1 



1 




with Qi = Q{li). Thus 



whence we find 



KL^ = — 

Qi 



Qi + Ki 

Ka Ql 



(6.22) 



(6.23) 



(6.24) 



It is now straightforward to take moments of (|4.54| ) to find {dK/dt)ioops- When we 
perform the r and ri integrations, we obtain a determinantal factor, which essentially cancels 
the normalization constants, and a factor corresponding to the expectation value of in 
the appropriate distribution. 

Let us consider the expectation value factor. In performing the r and ri integrations in 
the first term, it is simplest to change variables from r to r' = r + ri, so that becomes 
(r' — ri)^. Hence the appropriate expectation value comprises a sum of elements of the 
matrix L. In the second term, we have only a single element contributing. 

In this way we obtain 



'dK\ 



loops 



d/i A(/i 



I Ki,-K + 



1 rl 



Ql + Ki Jo Ql + Ki Jo Ql + Ki 



+ 



dyoKl 



From ( [4.40 ) and ( 6.25 ), we now obtain 



loops 



{ 2? V«'«'«, 
where I is another dimensionless function of the scale ratios defined by 

Qi{Ki-Aili)\ 



dlik{li) 2^/1 



Qi + Ki 



(6.25) 



(6.26) 



(6.27) 
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D. Rate of change of ( 



Estimating the rate of change of ( is less straightforward. The region of the integrand 
where li is very small is strongly suppressed by the A(Zi) factor, so we need not be particularly 
concerned with the behaviour of the remaining factors in this region. However, the same 
is not true of the region where the length / of our original segment becomes small. In the 
small-/ case there is a special feature that requires separate attention. 

We noted in Section [111 F| that for small values of li the conditional expectation value 
of is no longer given by the expression (|3.37|) but rather by (|3.51| ); it reduces to the 



unconditional expectation. The same argument applies to the conditional expectation value 
of for small /. It follows that in the contribution to (|6.25|) arising from the second 
[negative) term of ( |4.52| ), the integrand should in that limit be simply K ^ rather than 



K - K^JiQi + K 

The same argument applies to the first (positive) term. It is true that neither li nor / + li 
goes to zero in the small-/ limit, so the mean values of rj, r'^ or ri ■ r' individually should 
be well approximated by the Gaussian form. However, we are interested in the expectation 
value of r^, which is a small difference of these large quantities. Clearly, it too must approach 
/^ in the limit / — >^ 0. 

If we simply set the mean value of everywhere equal to K{1), the entire contribution 
would cancel out. However, this may not be quite correct. 

Suppose a loop is formed and that we choose a very small segment of length / on the 
loop. What is the variance of its extension? Clearly the leading term for small / is /^, but 
we could easily have 

[l + k)l^ 

K 



rLp = 1- . (6.28) 



with a value of A; 7^ 0. Indeed, it seems likely that loops are generally kinkier on small scales 
than long strings, in which case we would expect > 0. 
In this case, we easily find 

^ = k\. (6.29) 



VII. STRETCHING 

To complete our central task of deriving the various terms in the evolution equations for 
p[r(/)] and K{1), we have to estimate the various parameters and unknown functions that 
appear in them. For the stretching terms, we need to examine the parameters a and 
which are defined by (^191) and ( ^:20D . 

To match the notation of the preceding section, we denote the integration variable in these 
expressions by uq. Note that the expectation values appearing here are those conditional on 
the extension r of our chosen segment of length /. 
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A. Evaluation of [5 



First consider /5, given by ( [4.20D . To estimate the expectation value (q(yo)) conditional 
on the values of r, we again use the equation of motion for q, ( p.8|) . We could now derive an 
equation, similar to that for $ in the preceding section, but now for the complete angular 
distribution, $[q(|/i)] say, of q(|/i). 

However, we can use a simplifying Ansatz. The effect of r is important mainly in the 
region of relatively large values of yo — yi, in which it is reasonable to use a generalization 
of the exponential Ansatz (|5.14|) , namely 

where, as in ( ^.lOj ), 

b^3(q(2/i)). (7.2) 

Thus all we need is an equation for the expectation value (q(yi)). 
The required equation is a simple generalization of ( |5.27|) , namely 

(2 + Xy)^^^^ = H{piy,) - qq-p(yi)) 

+||(q(z/i)) - A(q(l/i)) + A(q(i/i))x, (7.3) 

with y = yo — yi as before. Here (q(i/i))x denotes the expectation value in the distribution 
X. Note that because p{yo) does not appear here, there is only a single stretching term on 
the right hand side; the factor of 2 that appeared in (|5.27|) is absent. 

Consider first the region of large y, where the expectation values are small. Within this 
region we may write, as in ( ^.31| ), 

(q(2/i))x = (p(2/i)) + Hvi)) (7.4) 
and replace (p(?/i) — qq-p(l/i)) by |(p(yi)). Then, as in ( p.33| ), we find 

(q(l/o)) ^ -^i^ fdy (l + 'e--y/'^'{piy,)), (7.5) 

where n is again given by ( 5.34 ). The factor | in ( 5.33|) is here replaced by | because there 
is only a single stretching term. 

The conditional expectation value (p(i/i)) which appears here was evaluated in the Gaus- 
sian approximation in Section HI. It is given by ( |3.41| ). 

As before, there will also be a non-linear contribution to (q(i/o)) arising from the sharp 
hole in the angular distribution for z close to 1. So long as the contributions to this quantity 
are both reasonably small (which appears to be the case), it should again be well represented 
by a sum 

(q(yo)) ~ (q(2/o))iin + (q(yo))ni- (7.6) 
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Moreover, the non-linear contribution to q(?/o)-p(?/o) is not much affected by the value of r, 
so it is reasonable to set 



(q(2/o))ni ~ -a„i(p(|/o)). 
Substituting from ( [7.61 ), ( fT^ ) and ([7.7|), we thus obtain 



(q(?/o)) = -anMvo)) - /f (l + ' ^-""'^"'^ 



\Piyo-y))- 



Now, by (g^ 



(P(.i)) = ^^jy^ r. 



Thus, (lOOl ) yields 



Ay' 



n-l 



We may now perform the integration over i/q, obtaining finally 

I" dy [K{1 -y) + K{1 + y) - 2K{y)] (l 



A + -H 

= ani + 



AK{1) Jo 



Ay 



n-l 



-xy/^e 



(7.7) 
(7.8) 

(7.9) 

(7.10) 
(7.11) 



(Recall that for negative values of the argument, K{y) is to be interpreted as K{\y\).) We 
note that at least in this approximation f3 is actually a function of / only, independent of r. 

Although we shall not require an explicit form, it is interesting to note that if K is 
assumed to have the two-scale exponential form (|3.10|) , the integral ( [7.11|) can be evaluated 
in terms of the incomplete gamma function. (Note that to do this, one must split the range 
of integration into separate ranges < y < I and I < y < oo.) 



B. Evaluation of a 

Now let us turn to the evaluation of a. The analogue of (|7.5|) is of course 

(p-q(l/o)) ^ l^^dy(^l + ^^ e->'y/'^\piyoypiyo~y)). (7.12) 

The factor of | is replaced here by | because, as in ( |5.33| ), we now have to include the effect 
of the rate of change of p{yo)- 

Substituting from ( p. 481 ) and adding a nonlinear contribution, we then get as before 

(p-q(yo)) = -am - f dy^^K^iy) + ^[K'iyo) + K\l - y,)] 

X [K'{y, -y) + K'il -y, + y)] ^ } + y) " ^ ^'""'"'^ ■ (^-IS) 
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First, let us examine the mean value a, given according to ( [4.19| ) by 

al = - f dyop^{yo). (7.14) 
Jo 

From ( |7.13| ), it is obvious that p^(yo) (to which the second term in the braces does not 
contribute) is in fact independent of z/q. Hence we find 

a = -Fq = a., + f dy K"(y) (l + y) " ' e"""'**' 



a„, + ^fdym (l + |^)""e-™«', (7.15) 



which may be compared with KC, equation (3.36). The important differences are the two 
extra factors in the integrand, especially the exponential representing the effect of intercom- 
muting, and the fact that A appears as well as H in the factor multiplying the integral. 
It will be useful to define the dimensionless function F by 

2.,„ = Fg|0 =A±i^f + (7.16) 

If / = ^K" is assumed to have the two-scale form (|0|) , and if the length scales of 
interest are large compared to 1/A, then F is expressible in terms of the incomplete gamma 
function: 

AH' 



with 



F^2w(^l + x-"e^T(n, x) , (7. 17) 



2B 

x=— + n. (7.18) 
A 



Returning to the remaining terms in a, we see from ( [4.19| ) and ( [7.13| ) that it may be 
written in the form (|4.16|) . The function a is given by 

= ^\tjn^ f dyo [K'iyo) + K'{1 - y,)] 



X 



K(l) Jo 

dy [K\y, -y) + K'{1 - yo + v)] U + y) e^^^/^«^ (7.19) 



Unfortunately, it is no longer possible to perform the j/o integration explicitly without 
assuming the form of K. 

C. Rates of change of length scales 

It is now easy to compute the rates of change of the various scale lengths. 
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The rate of change of L is given by ( [4.15| ) (together with ( [7.15| )). From (|6.15|) , we thus 
find 



sstr. TT /3 «nl 1 



For the other length scales, we examine the rate of change of given by ( [4.22| ). To do 
this, we need to examine the large- and small-Z behaviour of the parameters a and (3. 

First, we consider the limit I — * cxd, where K{1) k, 2^1 + constant. Substituting this form 
in the expression ( |7.11|) for /5, we find 



Pioo) = «ni + ^G, (7.21) 



where the dimensionless function G is given by 



A + |i/)frfWl + | e-A^^ (7.22) 



This function can be explicitly evaluated in terms of the incomplete gamma function, as 

G = 2 (^1 + n-"e"r(n, n). (7.23) 



If we use the same large-/ approximation in the expression ( |7.19 ) for la we see that the 
2/0 integral is of order ^^/. Clearly therefore 

a(oo) = 0. (7.24) 

These two results together yield 

^ = H[2P{oo) -a]=H (^ani + G - ^f) , (7.25) 

Next, we turn to the small-/ limit, using the approximation K{1) ^ /^ — l^/3(. Now in 
( [7. lip , the expression in square brackets in the integrand is clearly an even function of /, 
equal to PK"[y) + (9(/^). We shall need the terms in (3 up to order /. Since the denominator 
is of order /^, the /^ term in the integrand is irrelevant. The integral that appears here is 
then exactly the same as the one in ( [7.16[ ). 

It is no accident that /5(0) is almost exactly equal to a. In fact, it is physically obvious 
that for short segments, the extension r and the length / must expand by identical factors. 
There is an apparent difference between the two, namely the replacement of (A -|- |if) by 
(A + ^H). However, this difference is spurious. It arose because in the one case we included 
the effect of the rate of change of p{yo), while in the other we ignored the change of r. But 
for very short segments, r and p have of course essentially the same direction, and it is no 
longer reasonable to neglect r. Therefore in the small-/ limit, we ought to replace (A + |i7) 
by (A + |/7). 

In this way, we find 
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m ^ /So + = «ni + ^ + ^) (/ « C). (7.26) 

The / / C term arises from expanding the denominator. 

Next we turn to a. For small /, the leading term in the yo integral is 

' dy,[2l][lK"{y)] = 2l'K"{y). (7.27) 



Again, therefore, we find the integral F appearing: 

a{l) = ao + 0{l) = \F + 0{l), (/ ^ 0). (7.28) 

In this case, we shall not need the 0{l) term. 

To evaluate the term in ( |4.22| ), we need to examine the limiting behaviour of the 
function K(2){1)- As we saw in Section III F, as / — >• 0, K(2) oc /^/C; to be specific let us 
assume that 



K(2) ^ C-, 0. (7.29) 



Equivalently, C is given by 

K(2)il) CI 

1^0. (7.30) 



Now from (jO^ ), we find, using (fTW^ and {^^), 

= H{2>a - 2(3o + 6/3i - UCa^) = H(a^^ + F ] , (7.31) 



VIII. GRAVITATIONAL RADIATION 

There is one final effect that we have not so far considered, but which is in fact of great 
importance in the long-term evolution of the string network: gravitational back-reaction. 
We have been able to ignore it so far because it operates on a very different length scale 
from most of the other effects. 

Consider the gravitational radiation from a large length L of string. Essentially the only 
scale that can have any relevance here is the smallest scale C, which can roughly be identified 
with a mean inter-kink distance. We expect the rate of loss of energy, or equivalently length, 
to be 




gr.rad. 



where F is a constant of order 10^ or 10^. In other words, the lifetime of the small-scale 
structure would be of order (/TGfi. 
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It should be noted that numerically T here may be expected to differ somewhat from the 
values quoted in the literature, which have mostly been derived from studies of oscillating 
loops, because in those cases the length scale used was the length of the loop, whereas ( is 
defined somewhat differently. In fact, since the typical loop size is probably a few times (, 
our r is probably somewhat smaller. 

Another way of expressing is to think of the gravitational radiation as being gen- 
erated by each encounter between a pair of kinks, one left-moving and one right-moving. If 
we think of (, roughly speaking, as the mean inter- kink distance, we find that the number 
of such encounters on a left-moving segment of length L in a time interval dt is 

^. (8.2) 

Hence ( p. ID is equivalent to saying that each kink-kink encounter generates the release of an 
amount of energy equal to 

irG/i^C. (8.3) 

The gravitational radiation from infinite strings has been studied by several authors 
33| , |3^p7| . In particular, Hindmarsh has obtained a formula for the power emitted from 



encounters between left-moving and right-moving sequences of small-angle kinks. He finds 
(ref. ||2^, eq. (23)) that the power per unit length is 

r/P 

— = 2a2)G^^X0l Hd/TK)d'\ (8.4) 

where 9u and Or,, are the kink angles of the left- and right-moving kinks, d is the inter-kink 
distance and tk is the width of the string. Since the logarithm is slowly varying, it is 
reasonable to replace it with a constant, of order 10 to 10^. The length d is of course related 
to our C; by (13^31) , c/ ~ 

Now let us consider how to apply this formula to our problem, namely, how to estimate 
the rate of loss of energy or length from a chosen segment of left-moving string of length I 
and extension r. Consider first a relatively short segment, containing a single small-angle 
kink of angle 29. Clearly, we have 

- (8.5) 

Hence the formula (|8.4| ) suggests that the rate of loss of length should be proportional to 

As in the case of stretching there is an important consistency condition. The long string 
can be divided up conceptually into segments of any chosen length, and the proportional 
rate of loss of length must be independent of the choice, and must agree with (|8.1| ), i.e. 

1 r'^Ai; 

tgr.rad. — —I'- (8.6) 

Putting these two requirements together, we find that the expression for Zgr.rad. must be 
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''gr.rad. ^ T). T^n\' v"' ' / 



So far, we have concentrated on the change in the length / of our segment, but we should 
also ask whether there will be any change in the extension r. At first sight, it might seem 
that the answer should be no. Certainly in the case where I is large, gravitational radiation 
will change the extension at most by an insignificant amount. However, it is also clear 
that for small /, particularly if the segment is chosen to end near a kink, there could be 
a significant reduction in r. Indeed, as we shall see, it turns out that this is essential for 
consistency. 

If there is a reduction in r it seems reasonable to suppose that, for segments of a given 
length /, it too is proportional to (/^ — r^); certainly we would expect it to vanish in the 
extreme case of a straight segment, with = Let us therefore assume that it takes the 
form 



fgr.rad. = l^Hl)n ^771^^^ i^'^) 



where h{l) is an as yet unknown function. 

It also seems plausible to assume that for very large segments gravitational radiation has 
no significant effect on the overall extension. This would imply that 

h{l) ^ as / ^ oo. (8.9) 

We are now in a position to evaluate the change due to emission of gravitational radiation 
in the probability distribution p[r(/)]. It is obtained in exactly the same way as in the case 
of stretching, from the analogue of (|4.8|), namely 




= -^(/gr.rad.p) " ^ " (rgr.rad.p) + ^^^j^p. (8.10) 

gr.rad. 



From this, we easily find 



^ = -| (/gr.rad.r^) + 27.— Z + ^K. (8.11) 



Substituting from (|8.1| ), (|8.7] ) and (|8.8|) , and using the identity 



we obtain 



P-K P-K P-K' ^ 



As before, we can now find the rates of change of the various length scales. From (|6.15|) 
and (|8.1| ), we get 
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^gr.rad. _ TG/i 

Similarly, from ( [4.40|) and ( |8.13| ), we find 

^gr.rad. TGfl 



(8.14) 



(8.15) 



e C 

When we come to examine the third length scale, (, we can see why the presence of h is 
essential for consistency. Consider the limit I —>■ 0. For small /, the leading terms in K{1) 
are given by ( |3.50|) , so l"^ — K ^ ^^/3C- Moreover, as we saw earlier, 7^(2) ~ Cl^/(. Thus 
we find that the expression in the large braces in ( |8.13| ) has a term that behaves like for 
small /. This is inconsistent with the assumed form of K. Hence for consistency of our 
approximation we must assume that the coefficient of P vanishes, which requires that 



MO) 



1 



c 



(8.16) 



1-3C 

The value of h'{0) turns out also to be important, but is not constrained by any consistency 
requirement. 

Then, by (WM) , the term in (|8l^ ) yields 



gr.rad. 



c 



rG/i 



(8.17) 



C C 

where C is another constant, related in a somewhat complicated way to the leading terms 
in the power-series expansions of K, K{2) and h. Specifically, if 



K{1) = + 



/3 kX 



3C 12C^ 



+ 



HI) 



Cl^ 



c 



+ . 



(8.18) 



, / N hil 
h{0) - ^ + 



then 



C = -3(1 - 12C6 + ^hC) + 2/i(0)(l - QCe + 1:^0) + 6hi{l - 3C). 



(8.19) 



From (|3.56|) we see that if the model of Gaussian-distributed small-angle kinks is correct, 
then 



15 



(8.20) 



In any event, it seems likely that C -C 1. 

As a good first approximation, we may set C = and Cq = 0, which means that h{0) = 1 
and the value of /c4 becomes irrelevant. Then ( p.l7| ) becomes 



Cgr.rad. 



FG/i 



(8.21) 



c ^ ^ c 

Note, however, that the value of hi is clearly important. The sign of this term is crucial, as 
we shall see, in determining the nature of the solution. In our rough approximation, C > 
requires hi > ^. 
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IX. OVERALL EVOLUTION EQUATIONS 



We can now combine all the various terms together to give composite rate equations for 
L, for p[r(Z)], and for K{1). We hope to return to these equations at a later date. For the 
moment, however, we shall not write them down explicitly, but concentrate instead on the 
equations for the three length scales. 



A. Rates of change of length scales 

We begin with ^. From (fT^Uj ), and we obtain 



^ = H{ ---F] + -^E 



4 



1 



2^ 



(9.1) 



where the dimensionless functions F and E were defined in (|7.16|) and ( |6.17|) . 

For the other large length scale ^, we have, from ( |7.25| ), ( |4.41| ), ( |6.26| ) and ( ^.15| ), 



^ H (ani + G-^F 



we c ' 



(9.2) 



where / was defined in ( |6.27|) . 

Finally, we turn to the small-distance scale length (. Using ( [7. 311 ), ( ^4.431 ), ( |6.29| ) and 
( CT ), we find 



(9.3) 



where C was defined in (|8.19|). 



B. Estimation of parameters 

To determine the outcome of the evolutionary process, we first have to estimate how the 
unknown functions E,F,G,I depend on the ratios of length scales, and the likely magnitude 
of the additional parameters a^i, w, C, C and k. 

In some ways the most basic parameter is E, which governs the rate of loop formation A 
via the relation X = E/^. As we argued in Section VIB, if initially all the length scales are of 
comparable magnitude, then A will be very small because almost all loops formed reconnect. 
In that situation, E 1. Only when the upper cutoff ^/t has grown to be significantly 
larger than ( do many loops start to survive. Once that happens, we expect E to become 
of order unity. At present, we are not able to calculate the value very precisely, because it is 
strongly influenced by the intermediate-scale angular correlation effect discussed in Section 
VIB, for which we have only a very qualitative treatment. Fortunately, the precise value of 
E is not critical, because it appears as a common factor in several of the important terms. 

We note that, according to (|5.51|) , 
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This parameter will be very small, <^ 1, throughout most of the relevant region of 
parameter space, initially because E <^ 1 and later because ( <ti ^. There might be just a 
short period in which it becomes non-negligible, when ^ has grown large compared to ^, but 
^ and ( are still comparable in magnitude; but for the most part it can be safely neglected. 
Now let us turn to the function G which, in (|7.23|) , was expressed in terms of the incom- 



plete gamma function. It depends primarily on the value of A.^. According to ( |5.34 ), when 
A.^ is small, 1. Then we can use the fact that asymptotically r{n,n) ^ \^{n) to get 



2H\ 271 



TT 



^ ^ + 3A j V ~ V X^^ ^ '^^^ ^'"^^^ ^^-^^ 
In the opposite limit of large A^, where HC, is certainly negligible, we have 

G ^ -2Ei(-n) ^ -21nn ^ 41n(A0 (large AO- (9.6) 

Next, let us consider the function F, given by (|7.16| ). There are two large scales involved 
here (quite apart from the small scale (). The function F depends primarily on the two 
variables A^ and A^, though it also has a weak dependence on the small length scale. Note 
that, according to (|7.18| ), the argument x of the incomplete gamma function in ( |7.17| ) is 



2B 2w 2x 



Consider first the case where A^ ^ (-^O^- Then x ~ n and so (if H <^ X) 

F^wG, (Ae»(AO'). (9.8) 



Note that, directly from the integral definitions ( [7.16|) and ( [7.22|) , it follows that, at least 
when H is negligible, we always have F < G. 

In the opposite limit, where A^ ^ (-^0^5 have x ~ 2w/X^ ^ n, and consequently 
F <^ G. If A^ <^ 1, we can use the large-x form of the incomplete gamma function to get 

F ^ (A + \H)l ={l + ^) E (small AO- (9.9) 



3A 

On the other hand, if A^ is large, we can use the small-x approximation, obtaining 

F^2wln(Ae), ((AO'>Ae>l). (9.10) 

In the intermediate region, where A.^ and A^ are of order unity, F and G are both likely 
to be, very roughly, of order unity, with F < G. Strictly speaking, the expressions for F 
and G are both based on a linear approximation and cease to be valid for large A^ or A^. 
Recall that F/2 is the linear contribution to the value of a. From the simulations, we know 
that a ~ 0.2, so F is almost certainly significantly less than unity. 

Now we turn to the function / defined by ( |6.27| ). Initially, if all the length scales are 
comparable, /, like the other functions, will be small, because the loop formation probability 
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is small. Consider, however, a later time at which loop formation has started and ^ has grown 
to be ^ (. The integral is then dominated by the region where li is a few times (. In that 
region, Ki <^ 4^/i. Hence we have to a good approximation. 

Clearly, the ratio Q/K plays a very important role here. If it is approximately constant over 
the relevant range, then we get a very simple expression for /: 

Now, how large is Q likely to be? If the original expressions ( |4.59| ) and ( |4.60| ) for the 
loop-formation probability function 0(r, /) were correct, we should expect 

1 1 2a , , 

where the second term arises from the volume factor A. But in addition to this, we saw 
in Section VIB that the angular correlation effect would be expected to contribute another 
term to 1/Q. Although we have at present no means of estimating this contribution with any 
precision, it seems reasonable to expect that it too would be proportional to 1/P and of a 
similar order of magnitude. Even from ( |9.13 ), we see that 1/Q should be significantly larger 
than 1/K. The additional angular correlation term enhances this effect, suggesting that 
1/Q should be very substantially larger. Consequently, Q <^ K and it may be a reasonable 
first approximation to set Q/K = 0, which would mean that / ^ 2E. In any event, we 
expect that the ratio I/E will be not too far below 2. We shall see later that I/E is closely 
related to the parameter q introduced in CKA in fact I /E ^ q — 1. 

There are several other unknown parameters that enter our evolution equations. The 
parameter w is defined in terms of the large-/ behaviour of the variance function K: for large 
I, K ^ 2^1 — 2^'^/w. It was originally introduced in terms of the illustrative two-exponential 
model of Section IIIA, which suggests that it is limited to the range < w < 1; a typical 
value might be |. So long as it is not very small compared to 1, the precise value of w is 
not critical. 

The parameter C is defined by ( [7.29| ) in terms of the small-/ behaviour of the function 
K{2)(J')- Simple models of kinks suggest, e.g. as in ( p.20| ), that C is small compared to 
unity, say of order 0.1 or less. Within this general range, the precise value is probably not 
significant . 

The effect of gravitational radiation on the small length scale C, is governed, according 
to ( 8.17| ), by the parameter C*, defined by the rather complicated relation ( 8.19|) . In view 
of the number of independent parameters that contribute to this relation, it seems to be 
very difficult to estimate C from first principles. However, on physical grounds it seems 
clear that C should be positive, and presumably of order unity. The effect of gravitational 
back-reaction must be to smooth out the small-scale kinks on the string, and the expected 
time-scale for this process must be of order (/TGfi. 

Finally, we have to consider the parameter k. Like I/E, k is related to the parameter 
g — 1 defined in CKA |2T|, except that it is not concerned with the large-scale structure 
described by C, but with the small-scale kinkiness described by (. 
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One indication of its magnitude can be obtained by considering a related but slightly 
different parameter. Consider loops of length / and extension r. The probability of finding 
such a loop is proportional to p[r(/)]0(r, /). Hence if we choose a loop of length / at random, 
the variance of its extension is 

(9.14) 




K + Q 

Since we know that Q is small compared to K, this variance is <^ K. This statement applies 
to the whole loop, not to a segment on the loop, but of course the two are related. If a loop 
typically has an extension much less than that of a similar piece of long string, the same 
must be true for a small segment on the loop, though the effect is probably less dramatic. 
So we must expect r^^^^p < K{1)^ which implies > 0. It is not so easy to estimate its 
magnitude, but as we shall see there are reasons for believing that it cannot be large. 



C. Equations for scaling variables 



In our previous work to discuss the possibility of scaling, we expressed both our 
length scales as fractions of the horizon distance Rt. It turns out, however, that we can 
obtain slightly simpler equations if instead of the horizon we use the expansion time, 1/H. 
We define the three dimensionless ratios 

^'-''^ 

In the radiation era, 7 and 7 are identical to the variables used in KG and CKA, but in the 
matter-dominated era they are half as large. 

We also define p so that R oc t^/P, I.e., H = l/pt. Thus p = 2 in the radiation era and 
p = I in the matter era. 

The dimensionless functions F, G, E., I may now be regarded as functions of the ratios of 
these variables, e.g., 

i? = i^f-,-,-V (9.16) 

\e 7 p) 

Then, substituting into the evolution equations and dividing by if, we obtain 
7 /3 am F\ E VG^i 

-^^"-^+b-T-4j + 2^ + ^^' 

-pt^ = -p+ f ani + G-^\-^ + U + TG/ie, (9.17) 
7 V 2 / W7 2 

-pt! = -p+ (^ani + ^ 'l^^ ^) -^ + kE^ + TGfiCe. 

It is trivial to write down equations for the rates of change of the ratios of length scales. 
In particular, we find 

W7 



This will be useful later. 



54 



D. Scaling of ^ and ^ 



We now return to the main issue: what do the evolution equations tell us about how the 
length scales evolve? 

Let us suppose, as in Section VIB, that we start with all three length scales comparable 
in magnitude, and at least somewhat smaller than the horizon size. Initially, most loops 
that form will reconnect, so E will be small, as will the other dimensionless functions F, 
G and /. It is then clear that ^ will start to grow, because of the stretching effect, while 
initially ^ and ( will not. This will continue until loop production starts to be significant, 
when the upper cutoff, ^^/t exceeds a few times (. In that region, we have C ~ ^ ^• 

As A grows, it will first reach the point where ~ 1, while A^ is still small. In that 
region, E, F and / remain small, but G becomes of order unity. From ( |9.18| ) it is clear that 
the ratio C,/C will then start to grow. 

Some time later, we may reach the point where A.^ ~ 1 and A.^ ^ 1. If so, E, F and / 
would be of order unity but G would be even larger, so ^ would grow rapidly until it caught 
up with ^. 

Assuming the gravitational radiation terms are still negligible, do we then reach a regime 
where ^ and ^ at least approximately scale? 

For this to happen, the right hand sides of the first two equations in ( |9.17| ) must vanish. 
If we suppose for a moment that the values of E, F, G, I are known, then we can solve for 7 
and 7. The first equation always yields a positive value of 7: 

7 = i f2p-3 + «ni + ^V (9-19) 



E V "' 2 

The condition that the second equation yields a positive value of 7^/7 is 

{-2p + 2ani + 2G - F) + ^ (^2p - 3 + cini + > 0. (9.20) 

The ratio I/E plays the same role as g — 1 in CKA. Scaling requires a sufficiently large value 
of this parameter. 

If the condition (|9.2(]| ) is not satisfied, scaling cannot be achieved at the given values of 



E,F,G, L But of course these are not constants. If 7 satisfies ( |9.19| ) but ( |9.2U| ) is violated 



then clearly 7 will start to grow. The right hand side of the first of equations ( |9.17| ) then 



becomes positive, so 7 begins to fall. In other words, C, grows faster than C,. This of course 
affects the values of the various functions. U E = XC, remains of order unity, then A^ must 
grow, leading eventually to an increase in G (and a smaller increase in F). This in turn 
decreases the right hand side of the first equation in ( |9.17 ) and increases that of the second. 



leading us back towards scaling. It seems likely therefore that the parameters will adjust to 
reach the point where ^ and ^ do indeed scale. 

We must also consider what happens to the third length scale in this partial scaling 
regime. To find out, we examine the third of the equations ( |9.17| ). 



Let us assume that 7,7 <C e <C {TGfi) . Then both the intercommuting and gravita- 
tional radiation terms in the third equation are negligible. The essential question is: what 
is the sign of the right hand side? Using the partial-scaling value of 7, given by (|9.19| ), we 
find that for e to grow we need 
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k < ^ r-^. 9.21 



For example, if we take a^i = 0, C = 0, F = 0.4 and p = 2, then we require A; < |. (For 
a similar value of F, the condition is less restrictive in the matter-dominated era, where 
p = l) 

In passing, it is worth asking what would happen if the condition ( |9.21| ) were violated. 
In that case of course e would decrease; ( would start to catch up with C, and ^. If the 
dimensionless functions F,G,E,I were simply constants, C would continue to grow faster 
than ^ and until it actually exceeded them. This is obvious nonsense. What really happens 
is that eventually A starts to drop, thus returning us to a situation closer to our starting 
point. In fact, it appears that we would then reach a complete scaling regime even without 
invoking the gravitational radiation term. 

This should not be a surprising conclusion. We have already noted that k, like I/E, 
is analogous to the parameter g — 1 of our earlier work; k represents the fractional excess 
small-scale kinkiness of a loop as compared to a long string. If k is large, this means that 
the loop-formation process very efficiently removes small-scale kinkiness, so that complete 
scaling of all three length scales is indeed possible. 

What the simulations suggest, however, is a very different scenario, in which ^ and ^ at 
least approximately scale, but C does not. This is exactly what we expect if the inequality 
( p.21| ) is satisfied. 



Our conclusion here is really a straightforward generalization of our earlier results. In 



CKA we showed that scaling should occur if the parameter q — 1 exceeds a critical 
value. In our present work, we treat separately the large-scale and small-scale kinkiness. 
The parameter that represents the excess large-scale kinkiness is I/E. Provided it is big 
enough, ^ and ^ will scale (and, as we have argued, that is not a restrictive condition, 
because the parameters can adjust themselves until it is fulfilled). On the other hand, for 
( to scale as well, without becoming large, the excess small-scale kinkiness parameter k 
would have to exceed a critical value. This condition is almost certainly not fulfilled in the 
real system. 



E. Possibility of complete scaling 



Let us assume then that ( p.21| ) holds. Then while 7 and 7 settle down to scaling values, 
e continues to grow. 

This transient scaling regime continues until e grows to be of order (FG/i)^^, so that the 
gravitational radiation terms become important. The question is: do we then reach a new 
scaling regime in which all three length scales grow in proportion to t? Equivalently, do we 
reach a fixed point of the three equations ( p.lTI )': 



,? 



If there is a complete scaling solution, it is in the region where e ^ 7 ~ 7. This means 
that the intercommuting term in the third equation, XT^/^; is very small. If we neglect it, 
then it is straightforward, for given values of E, F, G, I, to solve the first and third equations 
for 7 and e, and then to find 7 from the second. 

The solutions for 7 and e are given by 
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{2p -3)C-p+{C + l)ani + + 3 - 12C)F 

(9.22) 

p-{2p- 3)k - (1 + fc)ani - 1(3 - + k)F 



e 



(C - A;)rG/i 



A consistent solution requires that the numerators of both these expressions be positive. 
This places a lower bound on C, namely 

^ )o-ani-|(l-4C)F 

— 2K ; ^ g_23 

2p-3 + a„i + iF ^ ^ 

where the second inequality is merely (|9.21| ). 

We recall that, according to (|8.171 ), C is the parameter that defines the efficiency with 
which gravitational radiation removes small-scale kinkiness. It is not at all surprising to find 
that this parameter must exceed some critical value if scaling is to be achieved. In fact, the 
second inequality in ( |9.23| ) ensures that loop formation is not sufficient to induce scaling of 
(, while the first inequality ensures that gravitational radiation is sufficient. 

The corresponding scaling value of 7 is found from the second of equations ( |9.17| ). As 
before, there is a consistency condition stemming from the requirement that 7^/7 be positive. 
This condition may be written 

{p-3 + 2ani + G){C-k)> (l-^){[2p-3 + «ni + lF]C - [p - «ni - |(1 - 4C)F]}. 



2E 



2* 

(9.24) 



Since both factors on the right hand side are definitely positive (the second in virtue of 
( p.23| )), this requires a minimum value of G: 



G' + 2a„i>3-p. (9.25) 

As in the earlier transient scaling regime, if this condition is not satisfied, then 7 will start 
to grow and 7 to fall until it is. 

We conclude that the parameters C and k are the most important in determining whether 
complete scaling will be achieved. The essential condition is ( |9.23|) . 



F. Stability 

Let us assume that a scaling solution exists and ask whether it is stable. To study this 
question, we consider small perturbations in the three scaling parameters, ^7, 57 and 6e, 
looking for solutions proportional to t^. Strictly speaking, of course, the functions E, F, G, I 
are functions of these parameters, but we shall assume that they vary slowly over the relevant 
range, so that they may be treated for the purpose of the stability analysis as constants. 

The condition for stability is that the three eigenvalues A are all negative. In other words, 
if we linearize the set of equations ( |9.17| ) about the presumed scaling solution, the resulting 



3x3 matrix must be positive definite. This matrix is 
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^ 

2X7 2Q_ _|_ L 

2X7 



2 



(9.26) 



There are three required conditions, the positivity of the trace, of the sum of 2 x 2 principal 
minor determinants and of the determinant. The first is automatic. The other two are 



—E + TGn 



C 



XT 



j -kE 


+ X 


7^ / 


>3 


+ -rG/i 


> 






g2 \^^^2 




e 





(9.27) 



and 



iC - k)TG^i— + ^ 

W7 e 



f4 + £_BU^ 

\W7^ 2 y W7e 



> 0. 



(9.28) 



In both cases, the terms are ordered according to the number of powers of the small quantities 
FG/i and 7/e that they contain. We see that in both equations, the leading terms are 
automatically positive. Thus we conclude that if the scaling solution exists it is almost 
certainly stable. The essential conditions are the bounds ( |9.23| ) on k and C required for the 
solution to exist. 



X. CONCLUSIONS 

We have analysed in considerable detail the various processes that affect the evolution of 
a network of cosmic strings — stretching, intercommuting, loop formation and gravitational 
radiation. We have derived a set of coupled evolution equations for the three length scales 
that describe the configuration — ^, the inter-string distance, S,, the large-scale persistence 
length along the string, and which characterizes the small-scale kinkiness. 

We showed that under reasonable assumptions and ^ reach a scaling regime, growing 
in proportion to the horizon, while C, grows slowly if at all. This continues until the ratio 
^/C becomes large enough for gravitational radiation effects to be significant. Thereafter, it 
is likely that a new scaling regime is reached, in which all three lengths scale, with ^ ~ ^ 
and C/e ~ TG/i. 

These conclusions depend on estimates of various functions and parameters that are 
still somewhat unreliable. The most significant parameters are fc, which describes the excess 
small-scale kinkiness of loops as compared with long pieces of string, and C, which determines 
the rate at which gravitational back-reaction smoothes the small-scale kinkiness. The fact 
that small-scale structure is seen to build up in the simulations strongly suggests that loop 
formation alone is not able to smooth the small-scale kinkiness, or in other words that k is 



less than the critical value, ( |9.21| ). The essential condition for a full scaling solution to be 
reached is that C exceeds the critical value, i.e., that gravitational radiation is effective in 
smoothing the small-scale kinks. 

It would obviously be desirable to be able to estimate the parameters k and C from first 
principles. We hope to do this at a later date. 

Another key part of the discussion concerns the build-up of angular correlations between 
the left- and right-movers due primarily to loop formation. We have given a fairly precise 
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account of the small-scale effect, leading to correlations between the p and q vectors, but 
the treatment of the intermediate-scale correlations involving segments of length a few times 
( is still rather sketchy. It is therefore very important to try to improve this part of the 
discussion. We hope to return to this question shortly. 

Once some of the remaining uncertainties have been resolved, we should be in a position 
to give a more accurate description of the scaling regime. It is very important to realise 
that the simulations performed so far, which neglect gravitational radiation, may have given 
misleading information about the typical scales involved. From the evolution equations, 
( ^.17| ), it is clear that when the gravitational radiation terms come into play, the effect 
will be to make 7 decrease. (It is easy to verify that the value of 7 in ( p.22| ) is less than 
( |9.19| ) provided that (|9.21| ) is satisfied.) In other words, the typical persistence length in the 
final scaling regime is probably a larger fraction of the horizon distance than it was in the 
temporary scaling regime where gravitational radiation was negligible — the only one so far 
accessible to simulations. Whether 7 also decreases is less obvious. That seems to depend 
primarily on the magnitude of the ratio I/E. We hope to provide a more detailed analysis 
of this behaviour in a future publication. 

There are many interesting questions that we should be able to address. With an analytic 
model it should be possible to estimate the typical density fluctuations, and hence microwave 
anisotropics, generated by a network of strings. Such an approach has previously been 
adopted in using the one-scale model, but it is essential that we understand the role 
the small-scale structure plays in determining the form and amplitude of the fluctuations. 
For example, do we see any evidence of the effects from the string tension renormalisation 
as indicated in the models of Albrecht and Stebbins 

Another interesting aspect which we hope to address concerns the transition period from 
radiation to matter domination. What happens to the scaling solutions in the two regimes; 
do they smoothly evolve from one scaling value to the other? A more detailed analysis of 
the evolution equations will enable a determination to be made of the approach to scaling. 
We have established the conditions required for scaling but have not discussed the time 
scales over which such conditions could be met. In particular, it would be interesting to 
examine the duration of the transient scaling regime and the way in which it is goes over 
to full scaling. This is important because the numerical simulations only directly give us 
information about evolution over a few expansion times. It might be possible to tackle this 
problem numerically in a simulation incorporating the effects of gravitational radiation. 

There are still many issues to be resolved, but we believe we have provided a firm 
foundation for future analytic work on the evolution of a network of cosmic strings. 
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FIGURES 
FIG. 1. Two contiguous segments 

FIG. 2. Three contiguous segments 

FIG. 3. Effect of intercommuting 

FIG. 4. Position of excised loop segment 

FIG. 5. Excision of a small loop 

FIG. 6. Coordinates on the string world sheet 
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